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APPLIED SCIENCES 


TWO-LASER RESONANCE ENHANCED PHOTOIONIZATION SPECTRUM OF GAS PHASE BENZENE 


Shanghai GUANGXUE XUEBAO [ACTA OPTICA SINICA] in Chinese Vol 7 No 1, Jan 87 
pp 1-10 


[English abstract of article by Fan Junying [5400 0193 4481] of Shanghai 
Institute of Optics and Fine Mechanics, Chinese Academy of Sciences; G. Muller, 
W.E. Schmid and K.L. Kompa of Max-Planck-Institute for Quantum Optics, Munich; 
J.L. Lyman of Los Alamos National Laboratory, New Mexico] 


[Text] This paper reports the results of a two-laser resonance enhanced photo- 
ionization experiment with benzene molecules. The excitation was provided by 
two frequency doubled dye lasers. The first laser pumped the molecule to a 
selected single state (S-1), from where it was ionized by a time delayed pulse 
of the second laser. The ion quantum yield depends on the intermediate vibronic 
state as well as on the wavelength of the ionizing laser. From the structure 
and intensity of the measured ion spectrum, the authors derived vibrational 
frequencies and molecular parameters of the ground electronic state of the ion 
that were highly precise. The contributions of autoionizing Rydberg levels to 
the ionization cross sections could be clearly distinguished from those of 
direct ionization. Some of these resonance peaks could be assigned to 
vibrations within these Rydberg states. 


The work was done at Max-Planck-Institute for Quantum Optics in West Germany. 
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APPLIED SCIENCES 


RESONANT SELF-FOCUSING OF TWO-FREQUENCY LASER BEAMS IN PLASMAS 


Shanghai GUANGXUE XUEBAO [ACTA OPTICA SINICA] in Chinese Vol 7 No 1, Jan 87 
pp 27-35 


[English abstract of article by Xu Tiefeng [1776 6993 1496], et al., of 
Shanghai Institute of Optics and Fine Mechanics, Chinese Academy of Sciences] 


[Text] In this paper the authors have studied theoretically the resonant 
self-focusing of two-frequency laser beams in plasmas. From the plasma fluid 
equations, the authors have derived an expression for the nonlinear dielectric 
constants of plasmas driven by two-frequency laser beams, showing the effect 

of the ponderomotive force of the plasma waves. After presenting an approximate 
solution for the longitudinal electrostatic field, characteristics of resonant 
self-focusing are shown. Finally, the influence of resonant sel f-focusing 

on particle acceleration in a beat wave accelerator is also discussed. 
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APPLIED SCIENCES 


STIMULATED ELECTRONIC RAMAN SCATTERING IN SODIUM VAPOR 


Shanghai GUANGXUE XUEBAO [ACTA OPTICA SINICA] in Chinese Vol 7 No 1, Jan 87 
pp 17-21 


[English abstract of article by Ma Zuguang [7456 4372 0342], et al., of the 
Laser Division, Harbin Institute of Technology] 


[Text] The first observation of the 3S-4S stimulated electronic Raman 
scattering (SERS) in sodium vapor with infrared output tunable from 2.38 ym 
to 2.65 ym is reported. The calculation for the power gain factors and the 
measurement for the infrared output as a function of dye laser tuning are 
made and are in good agreement. The pumping threshold as correlated to 
sodium vapor pressure is measured. The photo-conversion efficiency is about 
30 percent. In addition, the amplified spontaneous emission (ASE) of the 

Na 4S-3P transition is also observed. 
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APPLIED SCIENCES 


INTRACAVITY ABSORPTION SPECTROSCOPY WITH RING TRAVELING-WAVE DYE LASER 


Shanghai GUANCXUE XUEBAO [ACTA OPTICA SINICA] in Chinese Vol 7 No 1, Jan 87 
pp 11-16 





[English abstract of article by Yong Huayong [7167 5478 5554], et al., of the 
Department of Physics, Sichuan Universit: , Chengdu] 


[Text] By using a cw ring traveling-wave dye laser enhancement factor up to 
10° of the intracavity absorption, sensitivity was >Stained in air even if 

the pump power was two times over the threshold. sing the single-frequency 
scanning method for the first time, a resolution of 10~* angstrom was measured 
without sensitivity reduction. The experimental results agree with the 
predictions of the modified Brunner-Paul theory. 





























Fig. 2 Monitor and recorder for 


single-frequency swanning 

i—Dye laser, 2—Attenvator, J—X-Y recorder 

Fig. 1 Experimental set-up for broadband scanning (dashed-line é—1.5GHs scanning etalon, 5—C3-2100 oncillossope, 
represents the inserted parts in single-{requency scanning) 6—Mode] SMY generator 


1—Moade 360 Argon laser, 2—Dye jet, J— U nidirostiona! devios, 4—Thick etalon, 5—Thin 
etalon, 6—DBirefringent filter, 7—Galvanometer plate, 4—Lock circuit standard light 
resource, 10—Moiel (“T 50 monochrometer, 11— Photodiode, L2—X-Y recorder 
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APPLIED SCIENCES 


LASER INDUCED CHEMICAL ETCHING BY A SINGLE STEP 


Shanghai GUANGXUE XUEBAO [ACTA OPTICA SINICA] in Chinese Vol 7 No 1, Jan 87 
pp 49-54 


[English abstract of article by Li Ding [2621 0002], et al., of Shanghai 
Institute of Laser Technology] 


[Text] The scanning etching on Si, GaAs and Zn plates with a frequency 
doubling Ar* laser at 257.3 nm is reported. The width of the smallest etched 
line (1.2 m) has been obtained. The dependence of the etching depth on the 
light intensity/scanning rate is given. The diffusion coefficient, free path 
of the halogen radical and scale of laser focused spot have been estimated 
for the first time using the laser etching technique. A threshold phenomenon 
has been observed and explained. 
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APPLIED SCIENCES 





STUDY OF IR LASER PHOTOACOUSTIC SPECTRA OF ORGANIC MOLECULES ADSORBED ON 
METAL SURFACE 


Shanghai GUANGXUE XUEBAO [ACTA OPTICA SINICA] in Chinese Vol 6, No 12, 
Dec 86 pp 1057-1062 


[Article by Lu Huizong [7120 1520 1350], Chen Kaitai [7115 7030 3141}, 

He Maogi [0149 2021 4388], and Wang Zhaoyong [3769 0340 3057) of the 
Department of Physics, Fudan University, Shanghai; paper received 7 Dec 1985; 
revised manuscript received 24 Mar 1986; this paper was read at the 1985 
National Optics Annual Meeting] 


[Text] Abstract]: Usin: a branch-tuning CW CO» laser in the range of 

0.2 wm to 10.8 um we studied the IR photoacoustic spectra of organic mole- 
cules adsorbed on a silver surface. The absorbed molecular spectra of four 
layers of arachidic acid and cellulose diacetate with different surface 
densities was studied. No peak shift was found in a comparison between 

the IR photoacoustic spectra of solid arachidic acid near 944 cm™1 and the 
corresponding IR Fourier spectra of solid archidic acid. The IR photo- 
acoustic spectra of cellulose diacetate with o, = 1.4X10° cm™? and 

1 * 5.5x1015 cm=2 respectively was comr-red vith the corresponding 
transmission spectra of solid cellulose diacetate. It was found that the 
peak of the former near 1054 cm~! had a red shift of about 5 cm™! while the 
peak of the latter had no obvious shift within the range of accuracy of the 
experiment. 


I. Introduction 


Study of the physical properties of surface adsorption molecules has major 
significance for new chemical reaction catalysts, anti-corrosiveness of 
materials, surface oxidation properties of semiconductors as well as for 
material that have special behavior and surface treatment work. For 
surface research there are at present several very effective means for 
surface detection but they all has their limitations. The photoacoustic 
spectrum serves as a new detection means which can supplement the other 
methods to more clearly and effectively study surface physics problems. 
This is especially so for molecular adsurption. Because the vibratory 
state of characteristic molecular features generally are in ihe infrared 
region, infrared laser photoacoustic spectra will be an extremely effective 
means for the study of surface adsorption. Since the sensitivity of 











photoacoustic spectra is very high, it is entirely possible that it can be 
used to study monolayer and suhlayer adsorbed molecules.[1,2] This paper, 
using cellulose diacetate molecules as example, based on the adsorption 
coefficient of this sort of molecule in the veciaity of 1054 em71,13] pre- 
pared sample surfaces with densities of o,= 1.4x1045 cm72 2 celluloid 
series samples with specific gravity from 1.27~1.61 g cm™>, Experiments 
showed that photoacoustic spectral methods were quite sensitive when used 
to study the surface adsorption states. 


II. Experimental Apparatus and Measuremen 


The experimental apparatus was as shown in Figure 1. The light source was 
a homemade branch tunable CW CO» laser with wave length range of 9,2-10.8 um, 
maximum power of 10 W, and a linearly polarized beam, The photoacoustic 
detector was made from a S4 type piezoelectric ceramic transducer (PZT). 

One side of an optically polished slab (20 X 40 X 2) piezoelectric ceramic 
transducer was vapor plated in a vacuum (~10-6 Torr) with a high purity 
layer of silver about 1 um thick. Electrodes were led respectively from its 
two faces forming one detector. Because the sample on the surface is thin, 
during measurement after the silver base on the piezoelectric ceramic itself 
absorbs optical energy it produces a photoacoustic signal. This signal 
constitutes a background that cannot be ignored with respect to the entire 
measurement. 
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Fig. 1 Setup of the cxperiment 


1. Preparation of adsorption state sample 


In the experiment the preparation of the arachidic acid molecular sample 
was done with the die method to make an L-B film on the silver surface.|®] 
Arachidic acid molecules in a clean deionized water surface can form a very 
thin single layer film. If, along the direction of the film surface, there 
is applied a fixed surface tension, then it can form a dense monomolecular 
thin film.I7] The arachidic acid molecule forms a rod shaped molecule, one 
end of which is hydrophobic and the other hydrophilic. The hydrophobic end 
of the molecular film faces the outside of the solution surface. When 
preparing the sample, one slowly immerses the piezoelectric transducer into 
the solution surface and the hydrophobic ends of the arachidic acid 
molecules adsorb on the silver surface forming a monomolecular sample. Then, 
again slowly, removing transducer from the solution the first monomolecular 











film adsorbs another layer of which the hydrophobic ends face outward as 
shown in Fig. 2. By repeating this procedure, the metallic surface forms 
a polymolecular film. In these experiments four layer molecular samples 
were prepared. 





1 'ithin 





Fig. 2 
(a) Arachidic acid molecules adsorbed on the surface of evaporated silver film of the 
detector, which is being slowly thrusted in to the molecular film 
(5) Molecules adsorbed on the surface of evaporated silver film of the detector, which 
is being slowly lifted up from the molecular film 


To prepare the cellulose diacetate molecular samples we took 18 mg of 
cellulose diacetate sample solution and dissolved it in an 80 mg acetone 
solution. The, using a micro-dropper to takeup a sample, it was dropped 

on the rotating silver plating surface of the piezoelectric ceramic trans- 
ducer, Based on the molecular weight of this molecule and the surface area 
of the detector we estimated the molecular surface density of the sample on 
the silver surface. The surface densities of the samples prepared for the 
experiments were 0; = 1.4X1015 cm2 and 05 * 5.5x1019 cm=2, 


2. Method of measurement 


This paper adopted the differencing measurement method with a single light 
path, The laser beam was incident successively on two piezoelectric ceramic 
transducers as shown in Fig, 3. PZT(1) and PZT(2) are two piezoelectric 
ceramic transducers with the same properties, polished and silver plated 
under identical conditions then placed in a vacuum chamber as shown in 

Fig. 3. The surface of PZT(2) was prepared with the sample to be measured 
while no sample was prepared for the surface of PZT(1). With a p polarized 
laser beam incident on the sample of PZT(2) at an angle of 80°, the 
electrical signals produced by PZT(1) and PZT(2) were separately input to 

a dual channel phase locked amplifier. Doing subtractive measurements and 
simultaneously monitoring the power one can get 


, 


6(.) = SP2T) _ greta 


a 





in which gP2T(1) and gP2T(2) 


> are the respective photoacoustic signal strengths 
received by transducer £2). 





PZT (1) 


Fig. 3 Location of detectors. The angle 
between the two detectors is 160° 


Since the samples of adsorbed molecules are very thin films, the absorp- 
tion can be approximated linearly and considering that absorption of the 
metal itself is directly proportional to (1 - Rag)» in the end we can get 


__27a(d/c0s gp ) 
SO) 7 Bash, 9)" (1) 


where x is a constant related to the thermodynamic properties of the sample 
and the silver base, a is the absorption coefficient of the adsorbed sample, 
d is the thickness of the sample, and R,,(\,,¢) is the reflectivity of the 
silver base at wave length A, and incident angle ¢. From formula (1) we 
can see that S(\) is directly proportional to a(d/cos¢), that is in direct 
proportion to the sample absorption coefficient and also to the sample 
thickness d. Also S(\) is also related to Ry.(\,¢). However, from theoreti- 
cal estimates we know that in the vicinity of absorption peaks of the 
sample, the variation of (a d/sos¢) with wave length is much greater than 
the variation of [1-R,.(A,¢)] ~ with wave length. The reflectivity of 
metals in the tafeerell eaten and their conductivity at this time can be 
represented approximatley by their static values|5] 


Rug(A) 1-2/4, (2) 


11 fan 
[1-Ry(A)] 2V ce’ (3) 





A(i) = 





in which o is the conductivity of silver, A is the wave length, and c is 
the speed of light. For the wave length range measured in our experiments, 
9.2 um-110.8 wm, the corresponding (AA/AM™)) was 


4A/A(i) <0.08. (4) 


From formula (4) we can see that in the wave length range used in our 
experiments, since the variation of S(A) caused by changes in R,,(\,¢) was 
less than 8% while the absorption coefficient of general samples at 
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absorption peaks can be several times to several magnitudes more than at 
non-absorption peaks. The experimental results show that the influence of 
silver's basic reflectivity is very small and moreover that the change in 
this sort of effect with wave length is monotonic and therefore cannot have 
a very great influence on the spectral structure. 


In our experiments the modulation frequency of the laser was 23 liz and a 

small diaphragm (d = 2.0 mm) was used to obstruct a part of the edge light. 
The He-Ne laser in the experimental apparatus followed the same light path 
as the CO laser beam to act as an index beam for adjusting the light path. 


III, Experiment Results 


The experiment measured four iayers of arachidic acid, The results are as 
shown in Fig. 4, in which the vertical coordinate is arbitrary units repre- 
senting the strength of the photoacoustic signal. Figure 5 is the infrared 
Fourier spectra for solid sample of arachidic acid. A comparison of Fig. 4 
and Fig. 5 reveals that in the 924.90 cm! ~ 1086.84 cm™! (i1,e. 10.8 ~ 9.2 
um) band region, the peak at 944 cm7! in the photoacoustic spectrum corre- 
sponds to the absorption peak at the same location in the Fourler spectrum. 
The absorption peak of arachidic acid in the vicinity of 944 cm~! corresponds 
to the characteristic absorption of the carboxyl group cluster. Figure 6 
shows the measured results of a sample of cellulose diacetate with 
molecular surface densities of oj = 1.4x10) em7? and o7* 5. 5x1015 om”, 
Figure 7 is the infrared trnasmission spectrum of a sample of celiulose 
diacetate, The two spectra in Fig. 6 separately represent the photoacoustic 
spectra of the two samples with different molecular surface densities. From 
these results we see that the peak value amplitude and the surface density 
are in direct proportion which is in agreement with formula (1). 
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Fig. 4 IR laser photoacoustic spectrum of Fig.5 Infrared Fourier transmittance 
four layers of Arachidic acid molecules adso- spectrum of Arachidic acid 


rded on the surface of silver. The wavenumber 
of the peak is 944cem~! 
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Fig. 6 IR laser photoacoustic spectra of Fig. 7 18 transmittance spectrum of Cellulose 
Cellulose disectate molecules adsorbed on diacetate. The marked peak in the spectrum is 
the surface of silver with surface densities what we studied in the experiment. The wavena- 
of o3@1.4x 10" cm™ and o;=5.5x 1 mber of the peak is 1059¢m~* 
cm~*, respectively. The wavenumbers of 
the peaks are 1054cm™ and 1059cm™, 

respectively 


The peak value positions of the photoacoustic spectra are at 1054 cm -1 (co) ) 
and 1059 cm-l (c2) respectively. The former has about a 5 cm@! red shift: 
and the latter has no clearly observable shift within the experiment's 
range of precision relative to the peak position of this location in the 
absorption spectrum of Fig. 7. 


IV. Discussion 


That there is no observable shift of the peak value of arachidic acid's 
photoacoustic spectrum at the 944 cm! position with respect to the 
absorption spectrum is possibly related to the arrangement of arachidic 
acid on the silver surface, Arachidic acid is a long chain molecule with 
one end that is hydrophilic while the other is hydrophobic so that arrange- 
ment of its molecules on the silver surface is as shown in Fig. 8. In the 
figure the round head of the rod shaped molecule represents the hydrophilic 
end while the hydrophobic end,is in contact ~ the silver surface. The 
length of the molecule is 26 A and the 944 cm™+ peak corresponds to the 
carboxyl group culster of the hydrophilic end. ial Consequently the influence 
of the silver surface on the carboxyl group cluster is very small and the 
peak has no clear shift. On the other hand, the cellulose diacetate 
molecules cover the silver surface chaotically so there is always the ether 
and acetal structure of a portion of the molecules (corresponding to a 

1059 cm-l peak value) that will receive direct influence from the silver 
surface, When the sample is very thin this influence will be more pro- 
nounced since there is a relative increase in the molecules in direct 
contact yee the surface causing a red shift of the spectrum at the 

1059 cm * peak value. 


Photoacoustic spectra are a sensitive spectroscopy technique which can 
easily be applied in high vacuum conditions to study adsorption questions 
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of several simple molecules. The high vacuum environment can effectively 
suppress external acoustic interference but if one further applied low 
temperature cooling to the sample, then the effect of thermal noise would 
be reduced and in a high vacuum, low temperature environment we could expect 
a vast increase in sensitivity of measurement. Another course to increase 
measurement sensitivity is to use a surface electromagnetic wave method to 
make the laser beam and the surface adsorbed molecules more effectively 
interact. We used a prism coupled exciter surface electro-magnetic wave to 
excite the adsorbed molecules and then examined the photoacoustic signals 

of the molecule, i.e. surface electro-magnetic wave photoacoustic spectro- 
scopy. Preliminary results prove that the strength of photoacoustic signals 
produced with this method are several times greater than direct incidence. 
The detailed resuits of our work in this area will appear in another report. 
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Fig. 8 Arachidic acid molecules (two layers) 

stand up on the surface of silver, the sticks 

stand for the molecules and the dot heads of 

the sticks stand for the parts of positive 
hydrotropism 
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REFERENCES 

1, F. Trager, H. Coufal, PHYS. REV. LETT., 1982, 42, No. 23 (Dec), 1720. 

2, H. Coufal, F. Trager, SURFACE SCI., 1984, 145, No. 2/3 (Oct), L504, 


3. K. Bhasin, D, Bryan, J. CHEM. PHYS., 1976, 64, No. 12 (Jun), 15. 


4, “CRC Handbook of Chemistry and Physics", 53th ed., (CRC handbook 
18901 Cranwood parkway Cieveland, Ohio 44128 U.S.A.), C242, 1977. 


5. M. Bron, E. Wolf, "Principles of Optics", (Pergamon press 5th ed., 
1975, Oxford), 611. 


6, K. Boldgett, I, Langmuir, PHYS. REV., 1937, 51, June 1, 964, 


7, Georg Hass et al. ed., PHYSICS OF THIN FILMS, Vol 7, (Academic Press, 
1973), 311. 


8, Wang Zongming et al., SHIYONG HONGWAI GUANGPUXUE [PRACTICAL INFRARED 
SPECTROSCOPY] in Chinese, (Petrochemical Industry Press, 1978), 211. 


12966/9835 
CSO: 4008/34 13 











APPLIED SCIENCES 


GENERATION OF IR STIMULATED RADIATION, UV VISIBLE VIOLET COHERENT RADIATION 
BY TWO-PHOTON RESONANCE PUMPING IN K VAPOR 


Shanghai GUANGXUE XUEBAO [ACTA OPTICA SINICA] in Chinese Vol 6, No 12, 
Dec 86, pp 1063-1070 


[Article by Wang Zugeng [3769 4371 6342], Tang Xiaoling [0781 1420 3781], 
Zhang Kiachang [1728 7030 2490], and Zheng Yishan [6774 0001 0810] of the 
Department of Physics, East China Normal University; paper received 

4 Feb 1986; revised manuscript received 22 Apr 1986; first paragraph is 
source-supplied abstract} 


[Text] Abstract; When two-photon excited potassium atoms reach the 7S energy 
level, a series of photo-pumped stimulated emissions and cascade stimulated 
emissions are observed in the infrared region. The processes of photo- 

pumped stimulated emission and two-photon resonance four-wave mixing of 

pumped light produced corresponding strong coherent emission at 6P-4S and 
5P-4S. These two coherent emissions also had the function of pumped light 

and through a succession of four wave mixing processes produced several 
coherent emissions in the ultraviolet and the visible regions. This article 
discusses the process of generation of these observed infrared stimulated 
emissions as well as the ultraviolet and visible coherent emissions. 


Introduction 


The use of photo-puming to produce stimulated emission and new coherent 
emissions obtained from mixing has long aroused interest. In metal vapors, 
the two-photon excitation of Na, K, Ba, and Cs atoms through four-wave 
mixing and other processes has lead to observations of several coherent 
emissions, |1- 


In our work, we used a YAG laser pumped dye laser and with dual photon 
excited potassium atoms from the ground state, 4S, to 7S. An infrared 
detector received the photo-pumped infrared emission and a serics of cas- 
cade stimulated emissions from the 7S output. We used a photoelectric 
multiplier tube to detect these infrared stimulated emissions and the 

strong ultraviolet coherent emissions produced by pumped beam four-wave 
mixing interaction as well as to detect the ultraviolet and visible coherent 
emissions produced by these ultraviolet coherent emissions when acting as a 
pumping beam through successive four-wave mixing processes. 
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We identified the 10 infrared stimulated emission lines observed in the 
1.18 um to 7. 9 um region. We also discussed the formation process for the 
30-odd coherent emission lines in the 345.3 nm to 560.3 nm region. 


Experiment 


The experimental apparatus is as shown in Figure 1. A pulsed YAG frequency 
doubler laser, 1, pumps a dye laser, 2. When using DCM dye, the output 
energy in the vicinity of 660.0 nm is about 5 mJ and the laser output line 
width is about 0.01 nm. 








The dye laser beam is focused through a lens onto the center of the heat- 
pipe oven, 4, affixed to the potassium. The heat-pipe oven is 650 m long, 
with a diameter of 25 mm, and heating region lenggh of 250 mm. The incident 
and emergent ends of the pipe oven are fixed separately to quartz and CaF» 
windows, Using an auto thermostat the heat pipe oven is heated to 400°C 
where the corresponding concentration of potassium atoms is 1016 cm-3,[6] 


The directed infrared, visible, and ultraviolet emission signal emergent 
forward from the heat-pipe oven displays a light spot. At about 1 meter 
from the center of the heat-pipe oven, the radius of the light spot is 
about 3 mm, After passing through the filters, 5, the stimulated emissions 
and coherent emissions enters the monochromator, 6, for spectroscopic 
analysis. Appropriate band pass optical filters are used for the different 
signals in the infrared, visible, and ultraviolet regions to filter out the 
extraneous pumped light and leave the signals in the band to be analyzed. 
The infrared band signals emergent from the monochromator are received by 
a PbS detector of a 'Gaolai' pool while after being received by a photo- 
electric multiplying tube, the visible and ultraviolet band signals are 

fed through a signal processing system, 9, and recorded by the chart 
recorder, 10, and observed on an oscilloscope, 11. The detector, 7, 
received a minute quantity of pumped light to provide the trigger signal 
required by the signal processing system and the oscilloscope. 
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II. Results and Discussion 
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Fig.2 Energy level diagram 
of potassium and two-photon 
pumped stimulated radiation 


Figure 2 is a partial diagram of the potassium atom energy states.!7] When 
the output wave Jength, A; of the pumped beam is tuned to 660.45 nm, the 
potassium atom is excited by a two-photon resonance from the 4S ground 

state to the 7S state. The fixed direction emission spectral components of 
the forward emergent ray from the heat-pipe oven, after detailed measurement 
and computation can be divided into the following types based on the process 
of generation. 


1, Two-photon Resonance Pumped Stimulated Emission 


When the potassium atom is stimulated by two-photon resonance to the 7S 
state, population inversions form between the 7S state and 6P state and 
between the 7S state and the 5P state consequently producing 7.89 wm and 
1.80 um stimulated emissions, The stimulate emission coming from the 7S 
state also made the 6P state become largely population in turn forming 
population inversions between the 6P state and the 65S state and between the 
6P state and the 4D state thereby producing the 6.43 ym and 6.21 wm first 
degree cascade stimulated emissions. Similarly, second, third, and fourth 
degree cascade stimulated emission also is present. In the experiment, 
fourth degree cascade stimulated emission was observed as shown in 

Figure 2 and the final state of this emission was 4P. For greater clarity, 
the two-photon resonance pumping, the 7S to 6P stimulated emission, and 
the corresponding first degree of cascade stimulated emission processes can 
be represented as 


K (48) +20, —» K*(78), (1) 
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Ki(1S) —» K*(6P)+«(1S —6P, 7.89ym), (2) 


K*(6S) +o(6P -—6S, 6.43um), (3) 


K*(6P 
Kean) +w(6P—4D, 6.2ium), (4) 





Table 1 lays out the wave lengths of the 10 stimulated emissions and cascade 
stimulated emissions observed in the experiment as well as indicating the 
corresponding energy state transitions. 


Table 1 Stimulated and cascade stimulated radiation and corresponding 
transitions by two-photon, resonant pumping 









































No. corresponding trapstion or — 
1 18-062 7.89534 7.89 
2 ' 6P-+68 6.42252 6.42 
3 | 6P-24D | 6.21288 6.21 
4 4D-»5P 3.73412 8.73 
5 6S-»5P 3.66227 8.66 
6 5P-»8D 3.14154 8.14 
? | BF -958 2.70739 2.71 
8 18-»5P 1.90046 1.80 
. | 58 -»4P 1.25256 1.25 

10 3D -»4P 1.17729 1.18 











2. Coherent emission produced by two-photon resonance four wave mixing. 


From the above recounted four-wave mixing process by the interaction of 
two-photon resonance stimulated emission waves with the pumping waves, 
ultraviolet region coherent emission was produced. The form of mixing was 


wp = 20, — ca, (5) 


in which w, is the pumping frequency, We is the frequency of two-photon 
resonance stimulate emission or cascade stimulated emission, and wp is the 
coherent emission frequency produced by the process, Table 2 gives the 

9 coherent emission lines belonging to this sort of mixing that wer 
observed in the experiment. In the table, the first column is the coherent 
emission wave length values observed experimentally and the second column 
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is the values computed by formula (5). The third column gives the fre- 
quency mixing form, and the fourth column notes the relative intensities 

of each coherent emission. We note that iatensities were greatest for the 
coherent emissions produced by frequency mixing of two-photon stimulated 
emission waves and pumping waves (numbers 1 to 4 In Table 2.); with those 
produced by frequency mixing of first degree cascade emission waves and 
pumping waves next in intensity (numbers 5 and 6 in Table 2); while the 
weakest were those produced by the ‘“equency mixing of second (or third) 
degree waves (numbers 7 to 9 in Table 2); and the experiment did not ob- 
serve any coherent emissions involving fourth degree cascade emission waves. 


Table 2 UV colercnt radiation generated in two-photon resonant four-wave mizing 
































No. a ~~ form of mixing ralative intens.ty 

1 344.66 344.67 Quy — wy (7S — OP gn) 60.0 

2 54.76 544.74 Buz ws (7S —E6P yp) 76.0 

3 404.42 404.4) dg ws (75 — SPgn) 4.0 

. 44.76 4.72 Bury — wy (7S — BP yp) 61.0 

5 348.75 348.76 2ary — wg (OP — 4D) 18.0 

6 348.15 348.12 Mar, ~ ty (OP —68) 20.0 : 
7 863 00 32.96 Bary —ws (65 — SP) 6.5 

8 862.30 962.26 Qu; —w,; (4D—-SP) 6.¢ 

¥ 869.05 8¢8 .02 Boz — wg (SP — 3D) 4.5 

















From the polarization theory of nonlinear optics, the polarized intensity 
of the wp * Ju = We four-wave mixing process is 


Pop) = 80x (—wr, wr, Wt, — 9) EY (wx) B5(e), (6) 


Here Ey (uy) and Es (we) are the electric field intensity of the laser field 
and the stimulated emission field produced by two-photon excitation, and 
x 3) 4s the third order polarization frequency of the process 


x'(—wp, @1, @t, ~ wn) = K FS (7) 


Here N is the density of the atom, Ry, is the electric dipole matrix element 
between the ith and jth energy states, and Q, = [(Ei - E,)M) - ily is the 
complex frequency. [4 is a quantity related to the lifetime and relaxation 
of the ith energy state and K is a numerical factor. From formula (6) we 
see that the polarization intensity of the coherent field produced by four- 
wave mixing is in direct proportion to the polarization frequency x ) of 
this process, to the square of the laser field which participates in the 
frequency mixing, » and to the stimulated emission field, Eg. For the 
four coherent emissions marked as 1 to 4 in Table 2, there correspond re- 
spectively two resonance strengthening terms, i.e. 
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In addition, the stimulated emis~ion fields E,[wc(7S-6P)] and Es [we(7S-SP) } 
are the two strongest of the emissions in Table 1. Consequently these four 
coherent emission lines are very strong. The reason that emissions 5 and 6 

in Table 2 are second is because in the expressions for ,Q) we have 
respectively a single resonance factor, \ o~ 2 » and a nearly resonant factor, 
lgp-ep = 342 en) or Qep-wp = 292 cn] and because the first degree cascade 
stimulated emission intensities, Ecg[wo(7S-5P)]. The last two coherent lines 
in Table 2 (9 and 10) are very weak because in the expression for 

y{3), there is only a single resonant term and because the second and third 
degree cascade stimulated emission fields are even weaker. 


Figure 3 shows the two-photon resonance four-wave mixing form which 
produced the 344.6 m (6P 472-45) and 344.7 nm (6P) 72-45) coherent emission 


lines. 
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Fig. 3 Mixing scheme for generating 344 60m 
and 344.7 am coberent radiation Lines 


3. Second Degree Frequency Mixing Coherent Emission With wp as the Pumping 
Beam 


The strong coherent emission corresponding to the frequency wp * w (6P-4S) 
and op = w(5P-4S) produced by formula (5) can serve as a second degree 
pumping beam to produce four-wave mixing. At this time in the form of a 
single photon it pumps the potassium atoms and frequency mixes with the 
photon-pumped stimulated emission and the cascade stimulated emission 
described in the previous section. This produces a series of coherent 
emissions in the ultraviolet and visible region. The frequency mixing form 


is 


ig = Wp — Wy F Wy (8) 
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in which w' and wt are the stimulated emissions or cascade stimulated 
emissions listed Zn Table 1 and wu, is the coherent emission produced by 
the second degree four wave mixing. For example, ws = «.(6P-6S), 

we = ¥,(5S-4P), while up = 2u,; - #.(7S-6P), Then the four-wave 

alxine represented by formula (8) can produce coherent emission with wave 
length 513.5 mm. Figure 4 shows this sort of second degree four wave 
mixing. The polarization intensity of this process is given by the follow- 
ing formulas 


PP, — ex (— wo, wr, — wy — wt) Ep (wp) Ee wy) Be-(ws) (9) 


, . _N FpdioR alin 
x™(—wo, op, — &e, —w,)- XK > =? — opts —<éte . (10) 


For the frequency mixing form with a "+" in front of the we it is only 
necessary in formulae (9) and (10) to change the we to 
"we" and change Ec" (wi) to its complex conjugate Ecgu(we). 
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Fig. 4 Mixing scheme for genersting 
513.5 um coberent radiation line 


Table 3 gives the other 22 coherent emission lines observed in the experi- 
ment, There are produced by the four-wave mixing process with wp or uw, as 
the pumping beam. The third column in thetable gives the form of the four 
wave mixing which corresponds to the coherent emission observed. We can 
see clearly that all tzke w<(4S-6P) as pumping beam. Moreover, at least 
one first degree cascade st lated emission from Table 1 (1.¢€. we (6P-6S) 
or wo(6P-4D)) participates in the frequency mixing. So frequently not 
only are we able to observe in the experiment the corresponding second 
degree four-wave mixing coherent emission, but also the signal is a strong 


one (1 to 14 in Table 3). In fact, for conditions where wp(45-6P) secondarily 
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pumps the sample, the w¢(6P-6S) and w-(6P-4D) emissions also obtain 

resonance strengthening. In addition from formula (10) we know that for 
these 14 cohcrent emissions, their ,3) also have two resonance strengthening 
terms. The signals of the 8 ccherent emissions from 15 to 22 in Table 3 

are rather weak. From the frequency mixing forms to which they correspond 

we can see that some of them, although taking hang as the pumping 

beam, a strengthened first degree cascade emission did not take part in the 
frequency mixing, or that the two origina‘ resonance terms in the expression 
for ,9) for a "+" taken in front of wi in formula (8) now has become one 
resonaat term and one non-resonant term. Some take w,_(4S-5P) = 20;-wo(7S-5P) 
as pumping beam and obviously their emissions will be much weaker than with 
wp(4S-6P) = 2w, ~we (7S-6P). 


In summary, the above experimental facts are in very close agreement with 
theory. From formulae (9) and (10) we know that the strength of coherent 
emission fields produced in second degree four-wave mixing is in a direct 
porportion with pumping field strength and stimulated emission field strength 
of those that participate in the frequency mixing and with the polarizaticn 
frequency ~f the process. If the pumping or tic stimulated emission fields 
are weak the coherent field produced by the frequency will be weak. Also if 


Table 3 Coberent radiation generated through the second four-wave mixing processes 
(pumping beam: wp—2u, — we 79 —6P) or wp = 2, — wg 79 -- 5P)) 






































No. a ~~. form of mixing relative intens:ty 
1 345.30 345.27 @p — wy (6P —4D) +9 (6P-68) 12.0 
2 381.82 381.80 wp — wy (7S —5P) — wy (6P -63) 30.0 

a 3 382.60 382.57 wp — wy (78 —6P) —wy(6P -4D) 27.5 
‘ 386.10 386.09 wp — Quy (6P — 65) 25.0 

5 386. 90 386.87 wp — ty (OP — 68) —wy (6P - 4D) 28.5 

“ 387.70 set wp — Iw, (6P -41)) 25.0 

? 403 60 403 55 wy — wy (6P ~65) —wy (4D -5P) 25.0 

8 408 25 #05 29 @p — wy (OP —4D) — wy (65 — SP) 25.0 

9 412.00 411.46 @p — wy (6P — 68) - wy(SP - BP) 15.5 

10 412.90 412.86 tp — Wy (6P —4D) — wg (58 - 8D) 14.0 

FI 420.82 420.79 wp — wy (6P — 65) — ory (SP - 55) 13.5 




















[Table 3 continued on following page. } 
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[Table 3 continued, ] 





















































—_ 12 421.76 421.73 op~vs(6P 4D) —W@s(5P —5S8) 10.5 
13 513.55 513.54 wp — ws (6P —6S) —ws (58 —4P) 10.0 
16 514.90 514.87 wp—ws(5P-58) —wy(6P—4D) 8.0 
15 373.10 373.08 Wp —ws (6P— 4D) +H (SP —5S) 8.0 
16 481.15 481.20 wp —Ws(7S —5P) —ws(4D—5P) 4.5 
17 505.90 505 .93 wp —ws(78 —5P) —ws(5P —5S) 4.0 
18 531.62 531.58 wi, —ws (65 ~5P) —ws(5P~8D) 8.0 
19 546.45 546.49 wp —-Wg (6S —5P) —ws(58—4P) 8.5 
20 544.90 544.92 Wp — ws (58 —4P) +ws(“P—4D) 2.5 
21 560.30 560.27 w)p— wy (5P—5S) —w(5P—-3D) 8.0 
22 376.18 376.14 Wp — Ws (5S ~4P) +wg(78 —5P) 1.5 

SPI 
ws(55~4P) 


wp 


45 





6, (544 . 92nm) 








Fig. 5 Mixing scheme for 
generating 544.92 nm cohe- 
rent radiation line 


the resonance terms in the representation of the polarization frequency are 
few then the coherent emission field will be weak, 
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The reverse is also true. 














Figure 5 depicts the frequency mixing form producing a 544.9 nm coherent 
emission. Clearly, the coherent emission produced from this form of 


frequency mixing is weak. 


Figure 6 shows the particl coherent emission lines observed. The frequency 
mixing form corresponding to each spectral line is also indicated in Fig. 6. 


IV. Conclusion 


In this work we excited potassium atoms with two-photon resonance to the 7S 
energy state and observed 10 two-photon resonance pumped stimulated emissions 
and multiple cascade stimulated emissions in the infrared region. By the 
process of two-photon resonance four-wave mixing, 9 ultraviolet coherent 
emissions were produced, Of these the two strongest coherent emissions, 

wp = 2w; — wg(7S-6) and wp = 2w, - wo(7S-5SP), could also serve as a pumping 
beam to excite the potassium atom and by the process of second degree fre- 
quency mixing we detected 22 coherent emissions in the ultraviolet and visible 
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Fig. 6 Part of coberent radiation lines 
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[Text] Abstract; The paper reports a new Nd:YAG ultra short pulse laser 
system, The system's single pulse output energy at 1.06 um wavelength is 
greater than 100 mJ. Through a frequency doubler, pulse outputs of 50 mJ at 
0.53 um were obtained. The system was used on a range finder with main 
mirror of ¢600 mm at the Shanghai Observatory and measured the laser 
geodynamic satellite (LAGEOS). The ranging distance was 8000 km with an 
accuracy of 5 cm, 


I, Introduction 


Use of lasers for range finding of artificial satellites has developed from 
first and second generations and in the last few years, has advanced to a 
third generation. A major indication of this change in laser range-finding 
of satellites is the precision of measured distances, This precision depends 
on the laser pulse which acts as the emitted light source. The precision of 
the first generation was on the order of meters, the second was decimeters, 
while the third has reached centimeter magnitude. The U.S. Space Agency, 
NASA mobile stations all use several hundred ps locking laser third genera- 
tion laser satellite range finders. The lasers used are largely active 
locking lasers manufactured by the Xierfaniya company. The systems are 
large, expensive, and the switching components’ service life effects the 

life of the entire system. In the early eighties, pure passive, high 
efficiency locking lasers were already commercial products. Some stations, 
for example England's Greenwich Observatory and (Huoer) University 
cooperated|1] to use this system with a pulse width of 150 ps and single pulse 
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output of 30 mJ. Because a pure passive, locking pulse is constructed on 
the random noise of passive switching dyes in spontaneously radiated 
fluorescence, its locking probability is low, reproducibility is lacking, 
under repeated operation, the amplitude flutter between each time is large, 
around 150 percent, and often satellite pulses appear. 





In 1974 there was reported an intracavity added active modulator to improve 
the stability of passive locked lasers|8 adopting a high voltage 100 MHz 
signal modulated Pockels cell but its use is inconvenient. In 1978 an Nd 
glass laser with pulse series amplitude stability improved to +5 percent 
was reported, [3] Reference [4] reported a 1 Hz repetition operation Nd:YAG 
laser but did not provide detailed performance information, 


In 1981 we constructed a repetition active-passive Nd YAG locking laser 
using an acoustic-optical modulator as active locking component, improving 
the stability of the passive locked laser making the locking probability 
reach 100 percent and the amplitude stability at 10 Hz reach +2 percent. 

We did research on the cavity thermal stability for high repetition opera- 
tion, [5] making the repeatability of the oscillator reach 30 Hz. Through a 
single pulse selector, single pulse outputs of 0.5 to 1.0 mJ can ix 
obtained, and, through a two stage amplifier, at a wavelength of 1.06 um 
the output can reach 100 mJ. This ultra short pulsed laser system has 

been incorporated into the satellite range finder with main mirror aperture 
of ¢600 mm at the Shanghai Observatory as China's first third generation 
laser satellite range finder experimental system. Measurements with respect 
to laser geodynamic satellites (LAGEOS) show that its range extends to 

8000 km and its accuracy approaches 5 cm. 


II. Ultra Short Pulse Laser System 


co + 


9 
0 53,an 7 5 
“}<- - ey 
4 
Fig. 1 Schematic diagram of repetition rate ultrashort pulse laser system 
1—U)tre-short pulse oscillator, 2—Single pulse selector, 3, 4—45° mirrors, 


6,7, &—Beam expanders, 6—First stage of amplifier &—Second stage of 
amplifies 10—Harmonic geverator 























Our system is as shown in Fig. 1. The laser oscillator produces series of 
pulses from which the single pulse selector selects one to be passed through 
a two stage amplifier, The nonlinearity effects of a type II matching KDP 
crystal are used to do frequency doubling, putting out 0.53 um green light. 
Three doubling beam expanders are inserted into the system. The first two 
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of these are to expand the laser light spot to permit the laser beam to be 
amplified effectively as well as to lower the power density within the bean, 
preventing the nonlinear destruction of the laser medium at high power. The 
third beam expander is used to lower the power density to prevent damage to 
the frequency doubling crystal as well as to perfect the directionality and 
increase the doubling efficiency. 


The repetition ultra short pulse oscillator in the system has the features 
of high repeatability, high locking probability, and is highly stabile. 
Using the active locking set up, we performed studies of the thermal 
stability of the resonance cavity under locking conditions, [5] obtaining 

a set of computer modeled design methods that gave good agreement in 
experiments. 


The G factor of the resonance cavity is: 


Gs 


G3" Sots made ’ (1) 





in which a = a/b. a and b, respectively are the distances from the output 
lens and the rear cavity lens to the internal lens. With heating, the light 











spot dimensions on the front and rear of the laser ty and at the laser 
beam are 
W AL 2G? + 2aG, +a’ 
. x iig+a@ ° (2) 
AL 1 
W,-— - 
: x Gp+a’ G) 
W,-2k (2G? + 8aG, +a")? (4) 


@ (1+a’)(G,+a) (2G,+2aG,+e°) ’ 


model volume V —xlWj~ “iW [ (1-2) +(e) (5) 


Based on the major expressions above one can use a computer to make model 
computations and, based on the dimensions of the cavity's internal parts, 
select the normal operating parameters for the laser oscillator. 


Using an acoustic-optic modulator and dye respectively as active locking 

and passive locking components, active locking operation is attained. The 
medium of the acoustic-optic modulator is fused quartz 21 mm long with 

a transducer plate of Y36° cut LiNbO3. The transducer plate and the fused 
quartz are pressure bonded with indium. The dye is pentamethylene dissolved 
in ethylene dichloride solution. The dye box and the rear cavity lens 
comprise one unit which carries a circulation system to make the dye re- 
placement rate satisfy the repetition frequency. 


In order to lower the randomness and increase the reproducibility an active 


acoustic-optic modulator is inserted into the cavity to serve as a pre- 
selector. It forms an instantaneous window which has a width[6] 
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in which L is the cavity length, lj is the maximum single pass dissipation 
caused by sine modulation, x is the distance between the modulator and the 
nearest cavity lens, and N is the number of times the laser beam goes back 
and forth in the linear section of the cavity. From ~he above we know that 

tT gets small as x gets small, that is the modulator ought to be as close as 
possible to the resonance cavity lens. According to experimental measure- 
ments, x within a few cm, will not obviously increase the laser threshold 

and also will not significantly lower the stability. This is because the 
Nd:YAG laser is in a high gain state to within a few cm the constant dissi- 
pation action caused by x is minute. Formula (6) is appropriate to continuous 
locking lasers. In this paper formula (6) is used for a qualitative descrip- 
tion of the window width. 


The light path of the oscillator is shown in Fig. 2. In order to accomodate 
the echo time of the receiving system on the laser satellite range finder 
the pulse width is broadened from 30 ps to 120 ps using F-P to substitute for 
the common sphenoidal output lens and by adding a plate of auxiliary F-P 
changing the F-P thickness and reflectivity and altering the pulse width. 
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Fig. 2 High repetition rate ultra-short pulse oscillator Fig. 3 Mode-locked train of 


1—Rear mirror contacted by dye cell] with circulating system pulses (Time scale: 10 ns/cm) 
$—Transverse mode selecting aperture J—Xefiash lamp 
é—Froot mirror 


Figure 3 gives the envelope of a series of pulses. This is a 70 pulse re- 
peated exposure demonstrating that the stability of the pulse series is good, 


The pulse width was measured using a streak camera constructed by the Xi 

an Optical Institute. After the output was frequency doubled, it was split 
into two paths. One path was input to a high current tube producing a pulsed 
electrical signal to trigger the streak tube. The other path was input to the 
streak tube and immediately produced a streak signal, Figure 4 is the re- 
corded curve of the streak camera, The output light spot was a model TEMgo. 
Using a (nigro) meter to do a sweep trace with respect to the field diagram 
photograph, the sweep trace curve shown in Fig. 5 was obtained. 
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Fig.4 Pulse duration measured by streakcamera Fig.5 Scanning curve of output near feld spot 








We used an oscilloscope to measure a sample of pulse envelopes when the laser 
ocillator was operating at 10 Hz, revealing the pulse series envelope peak 
value stability. Figure 6 is a diagram obtained using a storage oscilloscope 
sampling 30 times at 10 Hz. From the 30 pulses in the figure the stability 
was computed to be about +2 percent. 





Fig. 6 Stability of pulse train (10 Hs) 


The output of the ultra short pulse oscillator is a series of pulses. To get 
a single pulse it is necessary to go through a single pulse selector (see 
Fig. 7). The light path is composed of two or more right angle prisms and 
one LiNbO3 Pockels cell. In the electrical circuit, the sequenced pulses 
are received by a PIN diode and after the produced signal is amplified it 
triggers a high voltage switching tube. A 2 ns leading high voltage pulse 
generated linearly by Blumlein shaping at a specially fixed instant is 
applied to the two electrodes of the Pockels cell. Of the sequenced pulses, 
the one that is jacketed by the high voltage pulse is the single pulse 
selected. 
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Fig.7 Single pulse selector 


1, $—Glan prism; 4—Pockels cell; é—Refiective mirror; 4£—Photodetector, 
6—Trigger-amplifer; 7—Forming network of high voltage pulse 


After beam expanding, the single pulse selected above is input to a two saage 
amplifier given a 1.06 um laser output up to 100 mJ. Through a type II KDP 
crystal frequency doubler, we can get a single pulse green laser of 50 mJ. 


The output energy of the single pulse of the Nd:YAG ultra short pulse laser 
system is adjustable 30 to 50 mJ (0.53 um). The pulse width is 120 ps; the 
type is TEMo90; with directionality 0.5 mrad, and repeating frequency of 

l Hz, 5 Hz, and 10 Hz. 


III, Application to Laser Satellite Range Finding 


The Nd:YAG ultra short pulse laser system combined with a main mirror $600 mm 
laser satellite range finder constitutes China's first third-generation laser 
satellite range finder experimental system. This system first measured 

laser geodynamic satellite distances on 12 Dec 1985. The measured data was 
sent to one of the international laser data processing collection centers--the 
U.S. space center, Goddard Space Flight Center Laser Tracking Net (NASA/GLTN). 
They sent the analyzed results back electronically. Table 1 is a summary of 
the actual observations of the ultra short pulse laser system and our analysis 
results. Figures 8 and 9 give the observed distance curve on 16 Dec 1985 

and Jan 8 1986 for the LAGEOS. On the entire curve there is a portion for 
which there is no data. This was the result of sudden malfunctions rather 
than any system defect. The analysis of Table 1 and 2 show that our use 

of the Nd:YAG ultra short pulse laser system brought the observation preci- 
sion to around 5 cm thereby attaining the level of foreign third generation 
range finding systems, 























Observation Date | Observation Time Length of are Observation Points <a 
Dec. 12. 1985 on~so— | iM 182 6.30m 
Dec. 16. 1985 20°12" ~41" 29 M 14 7.6cm 
Jan. 5. 1986 17%438@~56™ 13M 133 6.1¢0m 
Jaa. 8. 1986 17°27" ~ 378 10M 131 6.9°em 
Jan. 9. 1986 19916™~ 232 7M 16 5.0cm 























Table 2 NASA/GLTN Analysis Results 




































































Observation Date | COmrvation | Length of are | Observation | ciited Points [a 
Dec. 12. 1985 13%32~59= lM 80 19 72cm 
Dec. 16. 1985 20° 12"~41" 29 M 50 1 4.9¢em 
Jan. 5. 1986 1743" ~56= 13M 100 3 t.Scm 
Jan. 8. 1986 17°27°~37= 10 M 100 0 3.5cm 
Jan. 9. 1986 19°16"~23= 7M 16 2 6.0cm 
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Fig. 8 Range measuring curve of LAGEOS 
on Dec. 16. 1985 
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[Article by Cai Xijie [5591 1585 3381], Dai Meilan [2071 5019 5695], 

Lang Jiajun [6745 1367 0193], and Lu Ruixi [7120 3843 3556] of the Shanghai 
Institute of Fine Mechanics, Chinese Academy of Science and Qin Wenhua 

[6009 2429 7520] of the Department of Physics, Zhongshan University, 
Guangzhou; paper received 24 Feb 1986; revised manuscript received 7 Apr 1986; 
first two paragraphs are source-supplied abstract] 


[Text] Abstract: This paper gives the theory of tandem frequency doubling 
and relative orientation requirements of two crystals. Using two 1.4 cm 
thick type II KDP crystals and incident base frequency laser intensity of 
0.33 to 0.67 GW/cm2, the external conversion efficiency was over 60% and the 
laser beam diameter 42 mm. The theoretical results matched quite well with 
experimental values. 


We did comparisons of the experimental features of a tandem frequency doubler 
with L; = L2 = 1.4 cm to those of single crystal frequency doublers with 
thicknesses of 3 cm and 1.4 cm. The results showed that the tandem 
frequency doubler had the advantages of both the thin and the thick crystal 
frequency doublers. 


I. Introduction 


Due to many factors, in frequency doubling crystals there is always a 
certain phase mismatch present causing some difficulty in increasing the 
external conversion efficiency. Especially under conditions of low power 
density, only relying on increasing the crystal length to get the limit 

of the optimal length for a high conversion efficiency, as the crystal gets 
longer the reception half-width of the mismatch angle gets smaller so that 
correcting the phase mismatch becomes more arduous. People have adopted the 
idea of two crystals in tandem,[1~6,8] resolving this difficulty by making 
the mismatch angle the same in the two crystals and the signs opposite for 
the mismatch quantities which they cause. This can be achieved by changing 
the relative orientation of the two crystals. This paper gives the detailed 
theory for tandem frequency doubling. Using two 1.4 cm type II KDP crystals 
and a power density of 0.33 to 0.67 GW/cm2, we obtained a frequency doub|ling 
conversion efficiency of greater than 60 percent. 
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II. Theory 


As everyone knows the coupled wave equation describing the frequency 
doubling process is 


a+ 7:0, — —ik E°E,exp(—idkz), 


4, 
dz 


+5 7sB.~ —4K Bl exp (idk), and 


If we ignore the abscrption term, i.e. 1," ¥2° 0, and suppose that there 
is no dissipation inside the two crystals, then in a small signal approxi- 
mation, we can clea’ ly see the superiority of tandem frequency doubling. 


B,- xe" exp (3 [" 4b de ‘ae 


= ~4K Bi" exp (é4h2)de+-exp(é¢) |" exp (6 " ab de ar} 


1—erp(é i—exp(é 


in which ¢ is the phase difference between the fundamental and frequency 
doubled light in the crystal caused by the atmosphere and other media. 





When 4k2 = -Ak,, 1; = 1», 


By~ K BY 2S EEA) ) 1 + oxp (64)), (3) 
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EDP, KDP, 


Fig. 1 Relative orientations of tandem crystals 





Fig. 2 Principle of change in sign of 
phase mismatch 4k 
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If Akj = ok, it is only necessary for the relative orientation of the optical 
axes of the two crystals to be as shown in Figure 1. Figure 2 gives the 
principle for the sign change of Ak. For a single mismatched ray (like on 
the right side of the matching direction), because the signs of the mis- 
match angles caused in the two crystals placed as in Fig. 1, are opposite, 
since Sk = 846, thus the mismatch quantities, bk, and Sk, are oppositely 
signed. To obtain a very high conversion efficiency it is only necessary to 
adjust appropriately the separation between the two cyrstals to make ¢ = 2nr 
(the distance to produce ¢ = 7 in the atmosphere is 65 mm). At this time 


(4) 


Formula (4) shows that in the range of larger mismatch angles, the mismatch 
quantity of the first crystal is optimally offset by the second crystal. 

The result of the two crystals in tandem is this: the frequency doubling 
conversion efficiency is in direct proportion to the square of the crystals 
total length (L = 1) + 12 = 21,) while the reduction of conversion efficiency 
due to the phase mismatch Ak is just the same as for a single crystal 1). 
This lowers the requirements for adjusting precision and increases the 
frequency doubler stability. 


If ¢ = (2n +1)", from formula (3) E, = 0. 


That is, the difference frequency in the second crystal completely 
changes the doubled frequency obtained from the first ceyptal back to the 
fundamental frequency. This must be strenuousiy avoided. 4) 


In order better to understand the features of operation under high power 
densities and to approach as much as possible experimental conditions, we 
solved formula (1) numerically and did Gaussian numerical integration with 
respect to the time distribution of the light pulse. We also considered that 
after application of a refractive matching solution to the crystal surfaces 
there was still a 1 percent cross sectional refraction. The other parameters 
are found in reference [11]. Figures 3 and 4 give respectively the character- 
istic curves of incident light power density versus conversion efficiency 

for a single 1.4 cm crystal and two 1.4 tandem type II KDP crystal 
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Fig.3 Ds ibling characteristics of a 1.4cm type Il KDP 
crystal detuned 0~500 ured from phase matching; are 
raged over the Gaussian tempors! profile 
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Fig. 4 Doubling characteristics of a 1.4cm type-Ll KDP 
crystals detuned 0O~ 1000 rad from phase matching 


III, Experimental Apparatus 


We selected two 1.4 cm thick type II KDP crystals to use for the tandem 
frequency doubling experiments. The $70 mm KDP crystals were grown by the 
Chinese Academy of Science Fujian Structure of Matter Research Institute. 
In order to prevent deliquescence, the KDP crystals were sealed in 

a crystal cell unit. An antireflective film with respect to both the 
fundamental and doubled frequency was spread on the outer side of the glass 
window. A refractive matching solution was packed between the inner side and 
the crystal and to prevent destruction of the refractive matching solution 
at high power densities the entive KDP crystal cell units were assembled 
from ultra clean components. The total transmittance, T),, = 90 percent. 
The crystals were installed in two-dimensional servo adjustable frames with 


the angle adjustable to 1.5" per step. 
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The whole experiment was carried out at the 33 mm beam amplified output 

of the 6-path laser table in our laboratory. 9] The laser wave length was 
1.06 um; the pulse was 250 ps, the laser spot diameter was $42 mm, and with 
a plane wave front the directionality (measured using an array camera) was 
0.5 mrad in which the energy occupied 70 percent (0.6 GW/cm2 level). The 
light path of the experiment is shown in Fig. 5. We used (Ka) meters|7] 

1 and 2 to detect the input and output fundamental frequency light energies, 
El and Ej; (Ka) meter 3 to detect the output frequency doubler light energy, 
E2; filters 4 and 5, with Tj, = 70 percent, T,. = 0; filter 6, with 

Ty, = 9, Tz, = 68 percent; and beam eerresin, J 8, and 9 with R = 8 percent. 
The external conversion efficiency, Next = E2/F, and the internal conversion 


efficiency ";,, = E2/ (Ej + E.). 


i ‘ KDP KDP 
aw a ‘a ,5 49 
2 5 
6 


3 
Fig. 5 Experimental setup 








IV. Experiment Results and Discussion 


The debugging process for the tandem frequency doubler went as follows: 
starting with the first KDP crystal, it was tuned to the optimal phase 
matching direction. Then the second crystal was introduced and an angle 
tuning curve was made for the tandem frequency doubler with respect to the 
second crystal. After finding the optimal phase matching direction, a 
characteristic ".,+ - Ijg curve for the tandem frequency doubler was made. 


Figure 6 gives the tuning curve for the angle of the first crystal, with 
AS puny ~ 2400 urad, indicating that the mismatch half angle width A®pyyM 
is rather large. Figure 7 is the characteristic curve of the fundamental 
frequency power density versus the frequency doubler conversion efficiency 
for the first KDP crystal. In the figure, the theoretically computed 
curves for A® = 0, 600 yrad are also included. When the power density is 
1.33 GGW/cm2, we get the highest frequency doubler conversion efficiency 
in this experiment of Next = 45.5 percent corresponding to nj, = 50.5 
percent, 
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Fig. 6 Measured detuning curve of s 
1.4cm type-I1 KDP crystal 
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Fig. 7 Measured overall doubling energy efficiency of a 1.4cm-thick 
type-II KDP crystal as a function of input intensity J; with an 
sperture of 42mm. Dashed line: theoretical calculation for phase 
mismatch angle 49—(0.600 urad, averaged over Gaussian temporal profile 


Adding the second KDP crystal according to the requirements of tandem 
frequency doubling, Figure 8 gives the tuning curve of the angle for the 
tandem frequency doubler with this second KDP crystal. We measured 

A® = 2480 urad and the other values were the same as for a single 
crystal, This proves that in the matter of angle tuning, a tandem frequency 
doubler has the same advantage as a single crystal, that is A645), is fairly 
large. The required tuning precision is thus lower, 
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Fig. 9 Measured orerall doubling efficiency » of 
the tandem doubler (L,;=1,=1.4cm typeIl 
KDP) as a function of input intensity I, with 
an aperture of 42mm. Dashed line: theoretical 
calculation curve for phase mismatch angle 4j= 
0.600 wrad, averaged over Gaussian ternporal profile 


Figure 9 gives the characteristic curve of n for a tandem crystal 
frequency doubler. When Ij9 = 0.5 GW/cm2, we got the greatest frequency 
doubling conversion efficiency, n,,, = 63 percent corresponding to 

Nint = 77.8 percent, achieving a $25 mm wide diameter high efficiency 
frequency doubled output. In the range of incident power density, neither 
the refractive matching solution nor the KDP crystal itself suffered 
overload damage. 


In order to see even more clearly the features of a tandem frequency doubler 
we have collected together in Fig. 10 the experimental curves for the 1.4 
and 3.0 cm single crystal and the 1.4 cm tandem two crystal tandem frequency 
doublers. It is evident from this figure that the curve for the tandem 
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frequency doubler is better than that for the single crystal frequency 
doublers. At low pouer. ve < 0.3 GW/cm?), "tandem * "3cm? reaches a 
peak value at 0.3 GW/cm?2 le Ntandem reaches its peak at 0. 8° Gu /cm? - Both 
these results were ieasieiieadie predicted since the total length of the 
tandem frequency doubler, L = L] + Ly = 2.8 cm, is smaller than the L = 3 cm 
for a single crystal. The maximum conversion efficiencies ny ,,Gem and N30m 
were very close. However, it is noteworthy that in the higher power region 
(Ij9 > 0.5 GW/cm2), Ntandem ” . This is because the tandem frequency 
doubler has the capability of of tele ine the phase mismatch so that the fall 
in the conversion efficiency comes later for the tandem frequency doubler 
than for a single crystal of the same length. If we define the dynamic 
range of the conversion efficiency ng to be Q(no) = 11 0upperlimit/110lower- 


Limit | where the definition of li Oupperlimit and Ijo0lowerlimit 15 just for 


Ij in the region Tiqupperlimit 2 110 2 Tidlowerlimit» "7 No* From 
Figure 10 we know that 23¢m(55 percent) = 3 while Xr andem(55 percent) = 4 


so obviously the dynamic range of the tandem frequency doubler is greater 
than for the 3cm crystal. The actual performance of the tandem frequency 
doubler is good. By comparison with the 3cm single crystal experimental 
curve,[10] we can also see that the tandem frequency doubling, in addition 
to maintaining the rather high conversion efficiency advantage of a single 
thick crystal at low power, it also has a rather large mismatch angle recep- 
tion half width making the output of the frequency doubler more stable and 
the dispersion point smaller. 


Of special note is that the single 1.4cm crystal and two 1.4cm crystal 
tandem type II KDP freuqncy doublers both operate at An ~ 60 prad or even 
larger phase mismatch states but a single 3cm KDP crystal operates at a state 
of An ~ 270 uwrad so the difference between them is really quite large. 
However the tandem frequency doubler can still attain the same level of 
conversion efficiency, a point which completely explains why the phase mis- 
match offset capability of the tandem frequency doubler makes it require 
less than the single thick crystal frequency doubler with respect to laser 
performance, crystal material, and adjusting precision. A single 1.4cm 
KDP crystal, because it is thinner, at lower power density the conversion 
efficiency is not very high but at high power densities the conversion 
efficiency maintains a rising tendency, 
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Fig. 10 Experimental characteristics for 


three kinds of doublers 
o—L=1.4cm type-II EDP ($40402) , 49: 600 
urad~700 ured; b—L=3.0cm type-II EDP 
(840901), 49: 270 ured; e—L,—1,—1.4cm 
type-II EDP tandem crystals (840402), 4: 
600 ycrad~700 ured 


Finally, through comparison of these experimental data we reach the 
following conclusions: 


1. Single-thick crystals, although able at low power to attain high con- 
version efficiencies, at rising power density quickly exhibit a tendency for 
falling conversion efficiency. Moreover, since the conversion efficiency 

is more sensitive to the mismatch angle, the entire frequency doubler 

output is not stable. 


2. The conversion efficiency of single thin crystals, although not high 
at low power densities, in a wider power density range still maintain a 
rising tendency, are not sensitive to the mismatch angle, and thus have 
more stable frequency doubler outputs. 


3. Tandem frequency doublers combine the advantages of single thin crystals 
and single thick crystals. Not only do they attain high conversion effi- 
ciency at low power, they also maintain this high conversion efficiency with 
rising power density. Their frequency doubler output is also more stable. 


This work was undertaken with the good will and support of 

Professor Deng Ximing. The KDP crystals used in this work were provided 
by the Chinese Academy of Science Fujian Structure of Matter Research 
Institute. At that institute, Li Qingguo constructed the crystal box 
window multilayer film, Luo Yongshan designed the crystal box, and 

Fan Dianyuan provided the type number for the refractive matching solution. 
We extend our gratitude to all of these individuals. 


Also participating in the work were Shv Meizhong, Huang Kuixi, and 
Ge Xiaping. 
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GENERALIZED HYBRID VARIATIONAL PRINCIPLE, CORRESPONDING FINITE ELEMENT MODEL 
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Chinese Vol 7 No 5, May 86 pp 443-449 


[Article by Chen Wanji [7115 8001 0679] of the Dalian Institute of Technology; 
paper received 5 April 1985] 


[Text] Abstract: Based on the more generalized variational principle pro- 
posed by Professor Qian Weichang [6929 0251 7022], this paper provides an 
even more generalized hybrid variational principle suited for use in finite 
element methods. From this, a new general hybrid model is establishe”’. 


Further, with varying thickness thin plate bending elements as an example, a 
comparison is made of various hybrid elements as established based on differ- 
ent general variational principles. 


I. A More Generalized Extended Variational Principle in Elastic Mechanics 


Recently, Professor Qian Weichang proposed high order Lagrange multipliers 
and established a more generalized extended variational principle in elastic 
theory’ and consequently progressively unified and promoted study of a gener- 
alized variational principle. 


There are two places of insufficiency in the generalized variational princi- 
ple established in elastic mechanics using general Lagrange multipliers. The 
first is that one cannot from the potential energy principle directly obtain 
the Reissner principle. The second is that from the surplus-energy principle 
one cannot directly get the Hu Haichang [5170 3189 2490]-Jiujin [7750 3160] 
generalized variational principle. Evidently, the generalized variational 
principle set up using the former Lagrange multiplier method is limited. 


The high order Lagrange multiplier method proposed "xy Professor Weichang fur- 
ther developed Lagrange multiplier methods, overcoming the difficulties 
present in using traditional Lagrange multiplier methods to establish a 
variational principle. He also further established a still more general 
extended variational principle. 


The high order Lagrange multiplier terms proposed for Professor Weichang can 
also be introduced by penalty functions. This way is more direct and aug- 
ments the physical explanation for the high order multiplier terms. 
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When establishing a more generalized extended variational principle in elastic 
mechanics, the interior stress strain relationship is introduced as high order 
multiplier terms. In elastic mechanics, the stress strain relationship can be 
represented as: 


e=ac (1.1) 
or o=Ae (1.2) 
- ofe— ) e*Ac — L0ae=0 (1.3) 


in which €, ©O are,respectively,the stress and strain. a and A, respectively, 
are the flexibility and the elastic coefficient matrices and a™* = A. 

It is not difficult to prove 

rhe _ , ofa0 = : (e—ao)*(o—Ac)=— : (o—Azc)’a(o— Ac) =—(e—ac)*A(e—ac)/2. 
That is to say formula (1.3) is equivalent to using the stress grain coeffi- 
cient restrictions which are represented hy the high order terms. 





g7e— 


We take formula (1.3) to serve as a pnealty function and directly introduce it 
into a functional. Combined with theLagrange multiplier method we can derive 
two equivalent more general extended variational principles from the potential 
energy and surplus energy principles respectively. For the detailed deriva- 
tion see reference (1). Their functionals are 





eee ee 

-_ Tox=To.—({{A( ore e"Ac — L otae)uv eer 

in which Tomi peace ce— ore — Fru ]ar = [ftreds —frcu—aas (1.6) 
_ _ z ) _ _ 

Tay “{{{[ ore » the + (Do + )tu av — {er T)TudS — (1.7) 


r 5 


in which D is the differential operator corresponding to the equilibrium equa- 
tion, T is the boundary strength, u is the displacement, F is the given volume 
strength, T and u, respectively, are given boundary strength and boundary dis- 
placement on the boundaries S$ 5 and S,,, and \ is any given constant. 


Further, it can be proved that from 6llp, = 0 or dic, = 0, the variational 
equations corresponding to A taken as any constant are all equivalent to the 
elastic mechanics general basic equations. 


Obviously, Ty = “Ty 
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lip, and lic, are entirely unrestricted more general variational principles. 
Present general variational principles, including the Hu Haichang-Jiujin 
general variational principle, are special cases of this. 


Corresponding to different values of A, we can get various variational prin- 
ciples. 


(1) When A = 0, denote IIp) and IIc) as IIpo and IIe Obviously, this is 
equivalent to the Reissner p-inciple. Because Ilp) is directly derived from 
the surplus energy principle, these sorts of methods also are used to over- 
come another problem of the Lagrange multiplier method. 


(3) When A = 1/2, denote lp) and Ice, as pi, and Tigh, « Then we have 


Moy={{{{Ffrerae = +-0a0—o%e+207(D*u) |—Fru bay 


—{{truds—||T*(u—a)as 


Se 3. 
(1.8) 


Mey — HT ry 


This variational principle cannot be derived from the Hu Haichang-Jiujin vari- 
ational principle. If in formula (1.8) we substitute the restriction 
e = Dlu then 


IT py ={{f{ ; [ . e"Ac — + o7a0 + 0” (D*u) | - Prjor— []rress [frees 


; 5 (1.9) 


This is the general variational principle given by Liang Guoping [4731 0948 
1627) and Fu Zizhi [0102 1311 2535].’ 


If in formula (1.9) we substitute 


JN 0° (Dfu)dV ‘If. TudS— [ffcoeyrvar 


and also satisfy the restriction Do = 0, and do not count the volume strength 
ro 1 i = 1 —_ ? —_ —_ = 
IT py 2lJ{\( > e"Ac > ofac)dV {Jw (2T—T)dS {Js u)dS 


is js just the mixed variational principle used in reference (3). 


(1.10) 
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(4) When A = -l, denote lip, and IIc, as IIp-,and Iic-,. Then we have 


Ips {{ffetae—20re+.0" 0") +L oree—Fru]ay—[[truds—[ft%w—0045 Gs a1) 
7 Se 5. 


This variational principle also cannot be directly derived from the Hu Hai- 
chang-Jiujin general variational principle. If we substitute into formula 
(1.10) the restriction € = D'u (inthe interior), u = u (on the boundary S,) 


we get 


Ip..= j {ff (D?u) 7A(D*y) +L 0%ac—o*(0'u)—Fru Jer I! TTudS 


» 


(1.12) 


This is just the basic structure-potential energy principle given by Oden.“ 


When A is taken as different constants, one will get various new generalized 
variational principles thus making the various variational principles have a 
more unified framework. If we understand A to be the penalty function factor 
of the stress strain relationship restriction, then by the properties of 
penalty functions we know that when A + @ , the stress strain restriction 
will be satisfied. This way lip, and /-, will degenerate into the Reissner 
principle, the same as for 4 = 1. In approximation computations, the effect 
of different A values on the approximation precision is a subject awaiting 
further research. This paper will use an arithmetic example in the finite 
element method to discuss this question. 


Il. A More General Hybrid Extended Variational Principle in the Finite 
Element Method 


When IIp) and lic) are directly used in establishing a finite element model, in 
order to ensure convergence, commonly we must take into account the conformal 
requirements in Ilp, with respect to u or the conformal requirements in llc) 
with respect to 0. This is based on the fact that if Ilp) or [l-) use noncon- 
formal element functions, the establishment of finite element models is very 
inconvenient. If we relax the conformal restriction between elements for the 
Lagrange multipliers corresponding to Ilp, or !,,, then we can establish a 
generalized variational principle appropriate to the finite element method. 
According to the method of reference (5), we can establish two kinds of 
hybrid variational principles directly from Ilp) and lic). 


Hers= S| {|| fe7Ae—07 (ery) = Fruti(ote— eae - ; o*ao) |dv 
. 4p: 


- E JJ T'(u—a)dS—[(Trads (2.1) 
e Sf s, 
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or Maes = ff] — 3 erhe+ote+ (Do +F)tu—a (ore— Lethe 
. Ve 


_1 _ = . (2.2) 
2 cag ) av = JJ muse irut 


in which 0 is the displacement interpolation function on the element boundary 
3V,- It has two methods of construction. The first way is make & be uniquely 
determined from common parameters between elements. This is the hybrid method 
normally seen. The second way is to make U = (fu* + ou) /(a + 8), in which 

a, 8 are any constants (a + 8 = 0), with y® and u> being the displacement of 
the nonconformal element functions of the two sides of the element intersection 
boundary at the intersecting boundary. When a and § are taken with different 
values, i.e., a = 8 or a = 0, or 8 = O they are also frequently called parti- 
tioned general variational principles.*»’»* 


III. A More General Hybrid Model 


Based on the Il,p) or Il,yc, established models, we can call them a more general 
hybrid model. 


In the elements, suppose 


e=Ne 
o=PBthen T=RB (3.1) 
u=Bq 

G=Lq 





in which a, 8 are independent parameters in the elements, q is a nodal para- 
meter, N, P, B, L are the corresponding interpolation functions, and R is 
related to P. Substituting [l,.) in we get: 


To. L|- > a7Ha+6"Da +p" Wq—! (0TH 


—FH°Ch+pDa) — p°Gq | +470 


in which 


H= NTANGY (3.2) 


: 


4 


c=[{{prapay 


} 
D= Jere 
} 
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W= j |Jpirear 
G= | RPLGS 


vO {ff Prcqar + [frees 


From — =0 get <” —He + D°6 + A(Ha—D"S) =o (3.3) 

OTT acs : 3 4 

Fron ~—gg ~ =° get "Da + (W—G)q+/(C6—Da) =o (3.4) 

pron °u/=t* mo get (G—w)"B—0 (3.5) 

From (3.3) we get a=H-'D'g (3.6) 
Substitute (3.6) in (3.4) to get B= (1—2)DH-'D" + AC)" (G—W)q (3.7) 

LK*e=O 
Substitute (3.7) in (3.5) to get . (3.8) 


in which the element rigidity matrix K® is: 


K* = (G— W)*. (1; — 4) DH"'D* + AC)" (G— W) (3.9) 


When A = 0 we get the general hybrid element. When A = 1 we get the general 
hybrid stress element. As expected, this time the assumed stress satisfied 
the equilibrium equation in the elements giving the hybrid stress element. 


IV. A More General Hybrid Model for Various Thickness Bended Plates 


The hybrid elements relaxed the conformal requirement with respect to the 
element functions. Consequently, we can select a heuristic function from a 
broader range and make the performance of the elements be improved as much as 
possible. The general hybrid elements within the elements increased the 
number of variables and also relaxed as much as possible the restrictions on 
these variables. For example, the equilibrium condition within the elements, 
the continuity condition, the relationship of stress and strain, and the con- 
formal condition among the elements further improves the performance of the 
hybrid elements. 


For variable thickness plate elements, hybrid stress is inconvenient, making 
for the emergence of difficulties in the precise computation of the C matrix. 
Also seeking an inverse operation for the C matrix cannot be avoided. How- 
ever, the use of general hybrid elements can get od these problems.‘ 
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For the more general hybrid element, we give here a triangular thin plate 
element with linearly varying thickness and use a computational example to 
illustrate the effect of the A value on precision. 


The element nodal parameters were a single transverse displacement and two 
direction corners (Figure 1). 





\. @ 
| Sn 
. “n 





| a 


-_ ——_—-—- -—- —- =>» 





Figure 1 


The stress and strain within the element used linear interpolated values 
expressed by area coordinates. 


The displacement between elements used interpolation functions constructed 
from common nodal parameters between elements. The displacement w was a 
third degree function. The normal derivative dw/én was a linear function. 


From this we can derive the G, W, D matrix displays in the element rigidity 
matrix and H~’ C can be gotten through numerical integration. 


As a numerical example we computed for a cantilever plate with finite width 


and linearly varying thickness receiving a uniform load with dimensions as 
shown in Figure 2. The elastic constant was E = 2 x 10° kg/cm’, v= 0.3. 


a= ike om) | ff itil 
Din 
VARNERED A 


_——— "| 


LV 


= 106. 9a el” 
Af E 
VV" 


Figure 2. 
































Different values of A corresponded to different hybrid models. When 
A= 1, it was a hybrid stress element (the element stress field satisfied 
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equilibrium; when A = 0, it was a general hybrid element; when \ = 1/2, it 
was a hybrid element corresponding to the variational principle given in 
reference (2); when A = -l, it was the hybrid element corresponding to the 
variational principle given in reference (3). Im Figure 3 these are respec- 
tively denoted as R, H, L, and O. 








\ or . °Cr 
\ 
7 
wf 
at 
‘ j ——— 7 
Se -——_* 
cman SF SN 3 Is 
~— < 
: le i wien 
2.8. 2x20 
Se _— 
> “i ~-i @ € 2 . 
Figure 3 


For the computed results of the deflection for point A, see Table 1 and 

Figure 3. From Figure 3 we see that when the grid reaches a certain density, 
for different A values, the hybrid elements have nealry the same precision. 
For this type of element, if computational efficiency and saving of computer 
time is considered, the general hybrid element (H) is more easily implemented. 
From formula (3.8) we can see that for general hybrid element (A = 0) seeking 
the inversion operator can be avoided but for the other A values this step 
cannot be avoided. 


Table 1. Deflection of Point A (cm). Beam analytical solution uw, = 2.184 cm 








Grid 
<— eo, ix2 xa x6 x5 2x16 2x7 
A 
————e ———E = oo = = 
9 2.955 2.279 2.311 >.274 2.29 2.208 
i 2.728 2.8 2.320 2.271 2.26 2.206 
a | 2.637 2." 2.266 2.20 2.08 2.200 
0.5 2.000 2423 2.321 2.272 2.216 2.3 
~0.5 2.611 2.371 - 2.228 2.6 2.197 
2 | 2. 782 2.41 2.323 2.273 2.26 2.2068 
-? 2.986 2.388 80 2. 282 2.208 219 
5 2.990 2am 338 2. 2.223 2.20 
-§ | 5.00 2.264 2.20 2.211 2.200 2.1% 
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COMPACTNESS OF QUASI-CONFORMING ELEMENT SPACES, CONVERGENCE OF QUASI- 
CONFORMING ELEMENT METHOD 


Chongqing YINGYONG SHUXUE HE LIXUE [APPLIED MATHEMATICS AND MECHANICS] in 
Chinese Vol 7 No 5; May 86 pp 409-423 


[Article by Zhang Hongqing [1728 7703 1987] and Wang Ming [3769 7686] of the 
Dalian Institute of Technology, Department of Applied Mathematics; paper 
received 5 March 1985; first paragraph is source-supplied abstract] 


[Text] Abstract: This paper first discusses the compactness of quasi- 
conformal element space, generalizing the Rellich compactness theorem to 
quasi-conformal element space sequences with certain properties then it 
extends the generalized Poincare, Friedrichs, and Poincare-Friedrichs inequal- 
ity formulae to quasi-conformal element space. Then the convergence and error 
approximation of the quasi-conformal element method is discussed. The paper 
proves that if a quasi-conformal element space is approximate and strongly 
continuous, satisfies the element order condition and passes the IPT test, 
then the approximate solution is convergent. As examples, we show that con- 
vergence precision of 6, 9, 12, 15, 18, and 21 parameter quasi-conformal | 
elements under the L*>» 2(2) norm respectively are O(h_ ); Oth )> O(he ), O(he )» 
O(h? ), and O(h? ) order of magnitude. 


I. Introduction 


The mathematical foundations of the quasi~conformal element method were 
established in reference (1), where its convergence was summed up in three 
parts: 1) uniform positivity of the energy functional; 2) the approxima- 
bility of quasi-conformal element space; and 3) quasi-conformal element 
space passes the extended part interval test. Reference (1) discussed 2) 
and 3), but not 1). This problem will be discussed here. To reach this 
goal, we examine the compactness of quasi-conformal element space, general- 
izing the well-known Rellich compactness theorem to quasi-conformal element 
space sequences and extending Poincare, Friedrichs, and Poincare, Friedrichs 
inequality formulae to quasi-conformal element space. This way when utiliz- 
ing quasi-conformal elements to solve plate bending problems we can ensure 
uniform positivity of the energy functional. In addition, we improved the 
results of reference (1) concerning the approximation and interval parts 
test, proving that under rather broad conditions, the quasi~conformal element 
method converges as well as providing an error estimate. As example, 6, 9, 
12, 15, 18, and 21 parameter quasi-conformal plate elements converge. 
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As in reference (1), we take the fixed perimeter thin plate bending as an 
example to summarize our results. 


Let © be a bounded polygon in R*. Denote L*°7(9) = {u = (u°®, u?®, u°?, 
u**, w**, u°)| uiieL?(2), O€ i + j <2}, then L***(9) according to the norm 
below I 12.2(Q) is a Hilbert space: 


we€L***(Q), Iubis.eco = 5" [_ jaf? |Flx 
<a 
For \* w€H’ (2), define Ew= (Ww ,Ox, BW, Ox, u’, O,.2,u, O,0U. O5.2,%), then EH? (2) be- 
comes a closed subspace of L?97(9). On the space L***(2) define a bilinearity 


continuous functional a(+,*) as below 


u veL** (2). a(u,v) =|. (2797? + Day? yy?! 4 gy?) dx (1.1) 


For \/ f€(L?**(2))", consider the variation problem: 
w€H?(Q), of Fu, Fv) =f(Fv) (¥#v€H1(Q)) (1.2) 


Pick U,(t = 1, 2, ***) to be a finite dimension subspace of L*** (9), then the 
variation problem below 


u,€U,, a(u,.v) =f (v) (+#v7€U,) (1.3) 


is the multiple set function finite element approximation of (1.2). If U, is 
constructed according to the quasi-conformal element method, then (1.3) is 
the quasi-conformal element method to solve (1.2) (see reference (1)). 


For t= 1, 2, ***, let K, bea finite element dissection of { with the 
following properties: 


K, : for +K€ K,, K is a triangle (or a parallelogram), and Qe kK (r=t],2,---); 


Ko : the intersection or null OR of any two distinct elements in K, is a 
common edge or is a common apex; 


K; : denotes py to be the diameter of the maximum inscribed circle, hx is the 


outer diameter of K, hemsuphe (r=1,2,--,) then limA, =o and there exists a 


~- © 


positive number "1 >0 unrelated to T such that nh,<pe(¥KEK,.). 


Quasi-conformal elements construct U, according to the following 


u 
procedure,!»’»*» »° 
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First, given an integer t > 0, fory‘K€ K-,, given linear interpolation 
operators 17 }*.H*(K)-+H*(K) [._g: H'(K)+L*(@K), Wig: H'(K)+L* (aK) 204 175, 


H‘' (K)—~L*(@K) and the finite dimension subspace Ny J constructed from a 


multiple term formula, 0 < i+ j< 2. Next, define operator /7;', H'(K)~N;’, 


O<it+j<2. For WutHt (K), ye ju is determined by the following 
equations. 


+ PENY, [ PITgmdz=| | Pil.<wN ds —| a. P Meuds (1.4) 
¢ PENY, [PI wdx=| PIT, gwNN ds —| an PM gudx (1.5) 

+ PEN, [, PIT wdx= [PVE Sew —N NT ee) ds— | ,@ePIT tend 
| . (1,6) 

+PEN', 2{ PIywdx=| PANN JTS ews (Ni—Ni) M1 ew ids 

- | _(GayPIT fo +4x,P IT gw) dx (1.7) 

¢ PEN", [, Ply wdxe{ P(N; ew+N NT et) ds—| dnPIT'podx 
(1.8) 

in which N = (N,, N,)? is the unit external vector of 3K. 


Finally, define operator /7,, H*(Q)+L**(Q), for *w€H'(Q), [7.4 is determined 
by the following formula: 


(1 )) =I g(u\) (i+j<2, KEK.) (1.9) 


With respect to problem (1.2), make U,=/7,(H‘(Q)Hi(Q)), This is the quasi- 
conformal element method. 


According to reference (1), we say that (U.}, 41(Q) (orH'(Q)) possess 


approximability, if for ~weH!/(Q) (or H*(2)), we have 


lim inf [Ew—v,te22-9, 20 


¢-3 “ vw. Us 


We say that ({U,}, Hi(Q) (orl'(Q2)) passes the extended part interval test if 
when v-€U. (r1,2,-°). and supfe-tet20<™ , for @€C? (R*) tor'¥ pECG(Q)) 
we have jim 7,,,(g,0-)=0 (i@0,1,2,j—@ 1, 2). Here 


-o 
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ToalG ts) = [, (Gx uy +gul)dx 


Tosl¢ 2) =| (Ox.q 0% +") dx 
ra 


, 


Tialg,v.)= (Ox Ul + qui") dx 


- 


T12(9,%) = (xu + qu} dx 


- 


, 
Ts4(9.U-) = 9 Oueus tees dx 





Ts2(¢.0.)= | (Oapvt! + ut) dx 


With regard to the convergence of the solution of (1.3) we have the follow- 
ing theorem (1,6). 


Theorem 1. Let a(*,*) be uniformly positive definite, that is there exists 
a positive constant, @ unrelated to t such that 


YveEU,, ali rell*, , Sale, ve) (1.10) 


for tT = 1, 2, *** is uniformly established, then for 4/F€ (L***(9))', the 
necessary and sufficient condition for the solution u. of (1.3), under the 
meaning of |. | L*»*(Q) to converge on the solution of (1.2) Eu, is that {U, } 
and H§(2) is approximate and passes the extended part interval test. 


In the second section we prove the compactness theorem for conditions where 
{U,} possesses consistency and then prove the generalized Poincare inequality. 
In the third section for quasi-conformal element spaces with consistency, we 
prove that element order is able, under some conditions, to derive consistency 
and even stronger results. The major idea is the concept of normal affine 
continuous groups. The fourth section will discuss approximation and the 
extended part interval test, providing convergence theorems and precision 
estimates for the quasi-conformal element method. Last, section five will 
give proofs of the theorems in section three and section four. 


II. Compactness Theorems 


Let {U,} be a series of subspaces of L*°*(Q). We say that {U,} consistency 
if there exists an integer r > Pg! and a constant c unrelated to K and Tt such 


that: a) for4fv,E Ur, Yr, ViIIKEP,(K) for KEK, (05 i +4 < 2) is 
established; b) the rs Wt Me below 


le0,1, 20 vs \nese{ dy et lathe! Zot be} (2.1) 


teyel emrtel 
for\/Vr € Ur, KEK, \_ is uniformly established; c) the inequality below 


i+j<l, (03) (e) = (04) (x) |< ef A.sup a ler" (wy) | 
(®) mrartese (2.2) 


+hi*'"'sup sup > —_, , leew) 1} 


r ‘ 
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is uniformly established for ¥K’, K’€K,(x), +x€Q, +v,€l/, and +4¢r. Here 
(vi')*’. (vi')*" are the continuous extensions of vis|xr. vis | qe on K' and 
x”, K(x) = ik \x€K KEK.}, Kk’. K*. K are respectively the internal set of 
points of K', K", and K. 


We know that although U, is not necessarily contained in EH*(Q), it will be 
near EH*(Q) or its subspace. Therefore, the elements of U, must possess 
weaker properties similar to those of the elements of EH*(Q). Consistency 

is exactly what marks off these properties. The inequality (2.1) has in mind 
the relationship between the elements’ various components and inequality (2.2) 
pays attention to concinuity through interior boundaries. When U_ is not a 
quasi-conformal element space, the weak discontinuity in reference (7) and the 
weak continuity in reference (8) both are able to lead to consistency. The 
conditions for establishing consistency will be discussed below. Now we set 
up the compactness theorems. 


Theorem 2. If {U,} possesses consistency, then the following two conclusions 
are true: 


1) If v;€U, (t = 1, 2°+*) and under the meaning of L***(Q) it weakly con- 


verges to 0, then we have li . 
im v;’ 
"Poo toga! Ws" Tei =O, 


2) If {U,} and H*(Q) (or H§(Q)) possess approximability and pass the 
extended part interval test, then for any bounded sequence {v,} in L***(Q), 
v,€ Ur, there exists a subsequence {v; } of {v,} and vo€ H°(2) (or HG (2)) 

' os 
such that v,;'weakly converges to Ev» and in jo! —(E,)"4.4-0, = 0. 

ie ee oe bed | 

This theorem is the extension of the Rellich compactness theorem on the mul- 
tiple set function space U, and is also the generalization of the work in 
reference (8). Its proof entirely parallels the method in reference (8) so 
it is omitted here. Similarly, the methods of reference (8) can also give 
the theorems below. 


Theorem 3. If {U,} possesses consistency then the following conclusions are 
true: 


1) If {u;}and Hj(Q) possess approximability and pass the extended part 
interval test, then there exist constantsc, y such that the generalized 
Poincare-Friedrichs inequality 


Ce <e hs tv! —_ (440,€U,,) (2.3) 


is uniformly established for t > y. 


2) If {U,} and H*(Q) possess approximability and pass the extended part 
interval test, then there exist constants c, y such that the generalized 
Poincare inequality 
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CS <ef ) fu:’ lise t i ({,v: dx )} (2.4) 


tesud t+s<1 


and the extended Friedrichs inequality 


Weelis.tco) SL Do Wi Vag, +], lot? lds} (2.5) 


i+sad 


are uniformly established for 4v.€U,, ¥r>y- 


III. Element Order Conditions 


From here on we pick U, to be constructed by the quasi-conformal element 
method (see Section I). In this section we discuss the conditions for 

{U,} to possess consistency. We begin from the definition of element order 
conditions and a normal affine continuous group. 


Element order conditions. We say {U,} satisfies element order conditions if 
for YwEH'(K),4KEK, (r=1,2,---), from 7;°w=/7j'w=/7f*w=0 we are able to 
derive IT j°w [Tp wEP,(K), IT P°wEP,(K), and [T,cw—I[Th*wl ox =IlT}« w— 9 Tfwloe= 
G 

Mew gy Mi'w\se=0. 

The physical meaning of element order condition is "the zero energy mode is 
only able to be rigid body displacement" and the verification of this is the 
testing of the order of a matrix, so we call it the element order condition. 
Now we will discuss further the operators /7}*, /J,, , etc. For any element 
K, let there exist continuous functionals %.", ¢:,",°, ¢u.2, which are con- 
tinuous and unique on H®(K) as well as a set of polynomials Gye. Guys 
and a set of functions 9,,£, gi,£, 9%.e, 1<j<M, in the space {g|g€L'(8K), 
when F is a side of K, gir is a polynomial on F} such that for +w€H'(X), 


we have 
x 7 
IT?*w= 5 bs.a(w)G,,s 


j=! 


” 
IT, qgw= 3 by.0() O58 


jel 


= 3.1 
Tgqgw= 5 b5,0(w) o' 8 (3.1) 


gel 





wy 
Ty qwe 5 d,.2(w) 9", 


joi , 
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and when /7?*w=90_J],,w=/],,w=[15,e=0 B ¢;,4(w)=9,1 5M. We say *;,x is 


a parameter of pret 77, [3x, IT%x. 


We selected a fixed reference element K (a triangle or rectangle). For all 


elements, K, there exists an affine transform Fk = ByX + by, Bk isa2x2 
nonsingular matrix, beé R*, such that K = FyK. Thereupon for any function w 
defined on K (or 3K), denote 


(2)=w(F yt) (¥2€K(RAR)) (3.2) 


For the linear function 05 ,K» we define a linear function o; ,K on Ht (K) as 
follows: 


HEH *(R), b;.n(&) =d,,2(w) (3.3) 


in which i, W correspond to the form of formula (3.2). We say the operator 
order {/7%° Ion J7¢x.J7ixh=4{(1 8 Mex Mie, lox) |\KEK.,t=1,2,-} is a normal 
affine continuous group when: (1) if for a set of elements {X,};., there is 
a Bx, in the matrix norm sense converging to Bx Theu for j=1,---.M. ds.a5 
converges to $4 ,Kgunder the norm of the dual spaces of Ht (K), Cy Km uniformly 
converges to 654K,» while gj Ke? 84 3Km? and 8},K,> with the sense of L? (3K), con- 
verge respectively to 85,Ko° 855k, and aY x3 (2) for any element K and 
constant 8, denote R =m (5 |X Ox, ¥x€K}, and for any function w defined on 

K (or 9K), denote (5) =w(#-'x),xER (aR), then we have (/7$°w) (x) =(/7Pa) (Ox), 
(11yqw) (x) = (11g) (0x). (Ty ew) (x) =6 (17,2) (6x), (TF ew) (x) =6 (17, 5@) (0x), 

xEK (ROK), wEH' (kK), 


If the interpolation operators, /7?* J7,,, []jx./]%, only use functions of the 


interpolation points and their derivative values, and if certain kinds of 
integral average values of the functions and derivatives act as parameters, 
and if the relative positions of the interpolation points on K are affine 


invariant, then {/78*, [/yx,J7j¢, //3«! is a normal affine continuous group. 


This demonstrates that commonly used interpolation operators are normal 
affine continuous groups. It is not difficult to see that this is a 
generalization of affine groups.” 


We say that {Vz'}={N;'|KEA,,r=1.2,-"} is normal affine continuous if (1) 
for any element K there exists a group base P; x(l s & s Lys)» of niJ such 
that when Bx. > Bx,» we have Pa Kp converges uniformly to Pik, (l<sk<¢ 


L4;)s (2) for any element K and constant 6, if we denote K ee {5 |X Ox, ¥xEK}, 
then P,.,(x)=Pi,8(9%), (uxeK, 1<I<li;). 
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One special situation of {wid} being normal affine continuous is just that it 
be affine, that isthat there exists an wis such that for any element K, Nr = 


{P| P(x) =P (Fz'x). YPENYS. 

We say that {U,} possesses strong continuity if for eyery boundary F of any 
element K, there exist continuous functions qp, qe dp on the polynomial 

space which is defined on F. They have these properties ¢r(!)u;/')q7(1) ¥9,andigs} 
sige| OK KEK, ¥t}, {a3} @1ge|FCOK .KEK..¥4r} {gz} =iqiiFCOK, KEK,.¥,} 


respectively are affine groups such that when F=XK/)K’ is a side, for 
#wWEH'(Q), Gel Ta eto \ 9) =a9( LT oq 9), 95 (115 | 9) =@ 93 (i w\ 7), GF (13 ew| >) =-— 
(ve @\r), when F=KN@Q , for ¥wEH(Q) 1) H3(Q), ae(T.ew\r) =02( Mei!) = 
qi([1,4w\,)=0. We say {qr} is an affine group to indicate that if F = FyF, 
then for any polynomial P on F, we have g»(P?)=%¢(), and similarly for the 
others. 

If for ¥wEH‘(K), [lsew\r is used at the function interpolation points of 
win F, Tyee |r. and [T5 gw\» are respectively used at tangen- 
tial and normal derivative points when {U,} is strongly continuous; If 
[],~<w\,y is used at function interpolation points of w in F, and the integral 
average of /Ijqw\,, /1jgw\|, on F are uniquely determined by parameters of w 
on F, then {U_} is also strongly continuous, which is to say that the common 


spaces {u_} all are strongly continuous. 


Having the above concepts we can provide the conditions for establishing 
consistency. 


Theorem 4. Let the three following conditions be established: 


1) (172°, Qc. be. 13} is a normal affine continuous aw) and for 
¥YPEP AK), TPP=P, TegP wP\oe, The Pee ie, «P= oy lees 


2) (wy) is normal affine continuous and Ny contains a constant space 
(0< i+ j< 2); 

3) {U,} possesses strong continuity. 

Then from the establishment of element order we can assert that {U,) 


possesses consistency, moreover there exists a constant unrelated to F, K, 
and t such that the inequalities 


fm0.1, LO [My wl,nece LD Mews (3.4) 


teyed a+eet 








[85,17 {'w— 11 jw) ,,4+ |Pxu ll f'w— ITt'w\,.e <ch, X |Mitwl,s 


ivyet (3.5) 
[Myo Tege\.en<cht SL Mzwlys 
meeel (3.6) 
Ts gw— 2 18) +73 ww ° [7 2*w h . iM 
r as°°* ieee (4lge aN c owr™ec 7 = cle (3 7) 


for ¥w€H'(K), ¥K€K,, and x4, are uniformly established and when F=XK 1) XK’ 
is a side, the inequalities 


Te gw—1y_-@\,,.9< ht u | TF*wi.e + 77 2\ 40° ) (3.8) 


a+ te 


T3gw— 115, w\ 9+ | 3 pw— 1%. ~ . a 
’ -c oo? axe (Pleo rem cht > ( \l73*w\ e+ | 7, *\s0") 


(3.9) 


are established for ¥w€H'(Q) and when F=K/8@Q is a side, the inequalities 


| Tost ls.r< cht P7a ne (3.10) 
Diet leet Weel ir ch? 2 Mele (3.11) 


are established for *w€H'(Q2)1Hi(Q). 


Because these three conditions of this theorem are not indexed and commonly 
used elements satisfy these conditions, the element order condition can re- 
place the consistency condition. The verification of the element order con- 
dition can be summed up as an examination of a matrix order. 


For +w€H'(K), inN;’(i+j=2) select a group base, denote the coordinate of 
Il,’ in this group base as f}’(w) and make 

Bg (w) = ((Bp*(w))* , (By (w))*, (BE (w))*)* 

Oe (w) = (b5,0(W) 9,0 (Ww) baie (w))* 

Lelyt+lLy+Le (3.13) 


(3.12) 


in which Ly, is the dimension of ;'(), Then (3.6)~ (3.8) can be written as 
Agby(@) =Qrhx(w) (3.14) 


in which Ay is an L x L nonsingular matrix and Qy is an L x M matrix. 


59 





Theorem 5. Let ae NY contain a constant space and for ¥PEP(K), T}*P=P, 


TogP=P\og, 11:,P = rt on. 1% _ P= A lore Then the establishment of the 


element order condition is equivalent to the order of Qy being M - 3. 
IV. Conditions for Quasi-Conformal Element Convergence 


In the foregoing section we resolved the problm of consistency. Now we dis- 
cuss approximation and extended part interval convergence. First, the approx- 
imation result for affine groups (9) is extended to normal affine continuous 


groups. 
Theorem 6. Let ({J7f*. leg. Jig, 7%3e} be a normal affine continuous group and 


(niJ} be normal affine continuous, 0< i+j<«2. If there exists an integer 
ry 2 2 (1s i < 5), such that Dy Ho' (KR) -H*(K) (emer, +1), 

Nog Hts (K) LOK), ig Her (K) el OK) I1he, He '(K) LK) are continuous 
and for +PeP,,(K) we have J/f*P=P, for ¥PEP,.(K)we have TuP=P\o, for 
¥PEP,,(K)we have 7;, p - or te op. 1O° ¥PEP,,(K) we have 175, P= = lee, and 
Pry cae (KYCN (05+ 5<2), then there exists a constant c > 0 unrelated to K 


and + such that 


Oc lm, \w— IT ftwhieceket!~! lel oie (4.1) 
w= TT og? | 5. che 3 iwi, oye (4.2) 
Ow , r.—3 , 
as - Tyew| <ch,* Le (4.3) 
Ow r=} 
an —Miqw| - <ch, [le ene (4.4) 





are uniformly established for *w€H'(K),¥KEK,, t=1,2,-~ and 


S 
[Lu—T whim < eS he lel ss0 (4 


is uniformly established for Yw€H"(Q)0 H'(Q), ¥r=1.2.— . Here 


re max if,+l}. 

ise<e 
Because, under this theorem's conditions {U,)} is approximate and these condi- 
tions are rather broad, from now on we take. the condition that (U,} constructad 
using the quasi-conformal element method satisfies theorem 6 to eeen that it is 
approximate. 


Although under conditions of the energy functionals being positive and (U,) 
is approximate, passing the extended part interval test is a sufficient and 
necessary condition of convergence, its verification is not easy. Below we 
give an easier condition of verification--IPT testing. 








IPT Testing: For ¥K€K. there exists interpolation operators ITs... iT7, 
such that i//?. Neg. [1 %,. f,} are also a normal affine continuous group and 
the parameters of J] . Mae. ff:,. ff, are a linear combination of the para- 


. ° oP 
meters of MY. Mee. Wig, W%_ and for ¥PEP(K).z,P = 3?) MPa Oe 


If there exists r, > 2 such that the following two conditions are established: 


a) For ¥#€H'(Q)NHi(Q), when F=KOK’ is a side we have 


fT: u Tie 
¥PEP,,-2 (F), {, P| fs. jas=—[) | I |as (4.6) 
when FaK/@Q is a side we have 
¥PEPy-(F), | PMeuism| PMtpuds=o (4.7) 


b) For ¥wEH'(K), ¥KEK, (r=1.2,.--) we have 
Nt 15 ¢— 1% ,) @— NN, (113g — Mg) 
¥PEP y= 2(K), [P| 2N WN (1%, — M8.) + (NMI) (IT3g—Mtg)o | demo 
L Ni (115, = 9%, )w+ N,N, (11%, — Ete) , 


(4.8) 





then we say {U;} has passed the IPT test. 


Theorem 7. If {U,} possesses approximability and strong continuity, satisfies 
the element order condition and passes the IPT test, then for any /eé(L**(Q))’, 


the solution of (1.3) converges to the solution of (1.2): lim] Eu, —w, fos. ho =O. 


Theorem 8. Let {U,} possess approximability and strong continuity, satisfy 
the element order condition and pass the IPT test. If when F =K{)K’ is a 
side for ¥wEH'(Q), we have 


¥PEPn-a(F), | MoeuPdem| MoeruPas (4.9) 
established and when FeXK(\@Q for ¥wEH'(Q) 1 H1(Q) we have 
¢ PEP), -3(F), | Pll seudem0 (4.10) 


established then there exists a constant c, unrelated to u,, r such that when 
taking f(v) =| fords, fEL'(Q) , 1f mEH*(Q) H4(Q) then 


[Eu —welerrtca<e{ Soe tel ot o th? Itely } (4.11) 
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is uniformly established foryz, = 1, 2, **+. Here r' = max{r, +1, 4}, 
P_(F)=10}, rmmaxine+l.4h rf for KEA.. ¥t. Dew = DF weiss, ¥wEH"(K) 
u€H'(Q)9Hi(Q), Ul) Ee EC(®), viee=O). — thenr’ =ry+ 1.7 =maxirs+]). 


jaeee 


It is not difficult to verify that 6, 9, 12, 15, 18, and 21 parameter quasi- 
. conformal elements satisfy the conditions of theorem 7 (1-5). This proves. 
that they are convergent when used to solve plate bending problems. Table 1! 
provides their convergence precisions. 












Table 1. 
Wumber of |i ularity] yumber of : ty 
peremeters windies so Babee parameters _ ey aa 
~— _ --— - 
é Ole) ust 14*(0) 1S Om) wt HD) 
. O(m) ut H*(2) is O®) wt H*(D) 
12 0) ust H*(2) 21 ou) wt H*(Q) 

















Vv. Theorem Proofs 


Now we provide the proofs of the theorems in Sections III and IV. First we 
establish a lemma. 

Lemma 1. Let {279% . LJen, Dig. kg} be a normal affine continuous group and 
iw) be normal affine continuous. Then for a series of elements iAe«jié., and 
we€l/'(Ke) (m=g.j,---) if limi Ox. (Wu) ~ Gay (My) [=lim Bp0—Fa,1=°, we can 
assert that for 0 ¢i+j<¢ 2, (74 we) (Fg. #) uniformly converges to (/7}’ w,) 
(Faget). Here |+]) respectively are the Euclidean norm of the vector and the 
norm of the matrix to which it corresponds. 


Proof. Since ({//% , IJee, [1t,, I15,} is normal affine continuous, the conclu- 


sion of the lemma when i + j = 0 is obvious. For m= 0, 1, ***, in (1.4) make 
the coordinate conversion x = By_& + by, then use the differentiation chain 


rule and sum the following equalities 


dx = \det Bg, id#, dsm \det Bg,|{ B;' Vids, sary (5.1) 


and we can get: 
1< ll ely, \, Pie. Ty tin) (Fg &) (detBy \di 


=|), Pie. (Moratte) (Fut) (Bi! B),\detBe, ide 


7 (5.2) 
-|(o. 04, P ig. yay, tin) (Fg. #) \detBy, \d? 


e! 


in which 5;! = (67 ),,,. 
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; ia 
If denoting 1}, @= = LL FieP tas AL =| Pas. Py « \detBy,idz, R to be the 
right side term of (5.2), 1S 1, j © Lye, make 


Ag=( AY lie ha» Pu= (Pim) in XI! Ke= (i ) aes. } 


then (5.2) can become the following equality 
Am = = 


since {77¥ . [lus, ix. 11%.) is a normal affine continuous group and im") is 


normal affine continuous, therefore, by the hypothesis of the lemma, when 
mr 7A, * Ap, R, > Ro» and A, is nonsingular so 4, * 68). From this ve 


immediately derive that (/77, ~.) (F.#) uniformly converges to (J7¥w,) (Fast) . 
Similarly, other situations can be verified and the lemma is proved. 


Proof of Theorem 4: Proof of the longer parts is done in several steps. 


a) For ¥f>0, we prove that if the set J/,=(K|K is a triangle of the same 


area as K (or plane quadrilateral), 4«<€{}is not empty, then there exists a 
constant related to & such that 





fm0,1, S02 (MY whic ed) S [MFw),.. (5.3) 


eset at+eet 


[Om TP w — Tf @\,,.¢+ | Ox, 719 w~ 19 w\, e<elt) 5 [Tt*wl, (5.4) 
areci ° 
12 o— 1 oe0\ 08+ | 5-19 @ = Tigo) oe 
at sai | (5.5) 
+ | oc wm TE W | ee <c (€) > Vin le 


—_—- 
erect 


are uniformly established for ywEH'(K). ¥KE/;. We suppose that the above 
conclusions are not true. Then for m=1, 2, ***, there exists Ke€/; and 


w.€H'(K,), such that 


Lin Welite * > Tg. Wal te + |On iT? we~ TP. Daleks 


— 
eset 


+ 10x, 79 tam 19) We \etat | TE, tam Mine Mele0te 

+ (Mig, 5 1%, De let. + | 10g We =~", We ls.the 
(5.6) 
>m + iT We | eta 


seget 


63 











is established. Since for P€P,(Ka), 17%. P=P,[]* P=d:, P, iI. P=0;,P, 


—_ aP oP _ ~ 

te, P=S-|, » We, P= ayy, TH PHM, PHM, P=0, Mixa P=P lane, 
and these operators are all linear operators, for m=1, 2, *** we can let 
there be 


Ioxn(Wa)f=1, dxq (Wa) Lox, (/) (4+ PEP, (Ka)) (5.7) 
~ Simee {K atcJ., from reference (9) we know that {Bx is a bounded set. There- 
fore there exist series {¢¢,(wa)}, {Bxg,} (still denoted as {¢x_(Wa)},. {8x,, }) 

and element Ko and w,€H'(K,) satisfying. 


lim I9xn (Wa) —Ox, (%,) |= lim 18x.—Bx,l=0, Iox.(w,) | =1 (5.8) 


because By * BK» 055K, + $j,Ko. Thus for j = 1, ***, M we have limd san l)= 
" 2 2 
Osmo (1), lim $100 (0 (Beat): ) = b sine (x (Bx,%), )i=1,2) here By, = (bf )2.3 


it le! 


Consequently, from (5.7) we can get 


Ox, (W,) Ldxo(P) (4+ PEP,(K,)) (5.9) 


On the other hand, from the fact that ‘/7%, ox, [5,, 7%} is a normal affine 
continuous group and lemma 1, we know the left side of (5.6) is uniformly 
bounded. Thus Jim SJ" wa\.xn=0, that is (7", wa) (Fx, %) uniformly con- 


jeje 


verges to 0, i+ j= 2. Using lemma 1 again we get 
7% w= lI joss 
» “o U (t+ j=2) (5.10) 


From the element order condition we get //'* wo, 11) wEP.(Ko), 11% weEP (Ke), exo : 
— 0 

_ i) bw Welexe=0, Mix, wy 11% Wy | oxo 112g, Wy an J] Welexe=0, That is, if we 

make /, = im, Ww, then we have hy Wo =/7" P,, TT ox®y— 1] axe Pye IT, well ix, P= 

J1%,, we—113,, Py=0. Thereupon we get bxo(W,) =¢,,(P,). Again from (5.9) we get 

dy,(t,)=0, but this contradicts |¢x,(w,)J}=1 . Therefore, (5.3) (5.5) are 


established. 


1/2) Retk|R—0y x, x€K}. 


Then K has the same area as K and there exist constants c,, cz > 0 and 
unrelated to K, T such that 


b) For +A€K,, denote Ox = (area of K/area of K) 


cyhe* <Ox< che! , he<c, (5.11) 
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for +wEH'(K), denoting w(X%) = w(6;'x). +xER, Then since {J7%. Iss, J73,. 


17*,}, {NY } all are normal affine continuous, we have 





(117 w) (x) = Oy? (172 &) (Ox) (¥itj<2,x€K)  ) 

(17 exw) (x) = (17 o#@) (6 xx) (¥x€0K) (5.12) 
(175 ,w) (x) = Og (1755) (6x2) (¥x€0K) 
(113) (x) = 64 (I77,&) (64x) (+ x€0K ) 4 


are established and this way from (5.11) and(5.3)~ (5.5) we can derive (3.4)~ 
(3.7). Again, from (3.5) we know that (2.1) is established. 


c) Now we prove that for £ >1, there exists a constant c(&) which is 
related only to € such that the following inequality 


Dox w— Tox\yr+ | Tigw— Migr w\ye + (ITS, w— 5p wy (5.13) 


<e(é) S {UT wie + 1D) whee} 


t+jae 


is uniformly established for ¥w€//'(A UX’) and this sort of K and K': 
F=K1\K’ 4s a side, K and K are the same area, fig<t, hg<t, py >t hy’. 


If (5.13) is true then using a method similar to b) we can prove (3.8)~ (3.9). 
Now suppose (5.13) is not established, then for m= 1, 2, *** there exists 


kK KD possessing the properties K,,{/|K%=F,, is a side, kK and K are the same 


area, hy <t, he’ <i, ps’ Do thy! Mk weEH' (Ke KL) such that 


| Teng Wa — Loy’ Wa\ s.r. + Tix, Wa 14, We lo teat | sy wn 


<5 2 alte > So CIMY, alte + |TY; waleste 3} (5.14) 


t+jed 
and, like a), we suppose for m= 1, 2, ***, we have 


idee (Wm xe)!" + lox: (We, ige b=), bee (Waiz,) Lbce(P ie.), (5.15) 
¥ PEP, (Ka J K2) 
and there exist Ko, K'o, and weéH'(K,\/KA4). F.=K,9K{ is a side, such that 


when m*> © 


7 a Ke) ~* ao (ty ae! ys { ie i he )->dx (wyixy) 


bs. Bay, Bus -» Br! | (5.16) 


and we have 


1 w,=0,/1") w,=0 i+ j=? 
Ped , (i+j ) (5.17) 
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Thus, we have Gne(Meixe) = 0, 1% wEP (KS). Mex. Dy— LI. 9, lone =lI% 5. > — 

C) 
ds’ 
qF, (17351, Wel re) = 9? (1%, We |r) =0, thus, {T*: wEP,(KS); and Gr lllex’s Weise) = 0, 
so iT. w,=0. That is x’, (w,) =0, which is in contradiction with 


, 0 
i]*: We | ex; =I] fe; Wy — saz, 11". Weieg, =O. From strong continuity we know 





idxe (Wei xe) 1*+ idx, (wla)f*=1 ° Therefore, (5.13) is established whence (3.8) 
and (3.9) are established. 


Using (3.5) ~(3.9) we can prove that when F=K1K’ is a side, the inequality 


aod | Tew — 11. wie the? TY w 
t+jei 


- (5.18) 
— fy", wisgeQe ) CLT oe + | IT? w| 0°] 


is uniformly established for ywEéH'(Q). Here c is unrelated to K, K' and 
K, K'€ K,. Using (5.18) and the lemma of reference (1) with the method of 
reference (8), we can prove formula (2.2). 


d) Last, using a method similar to that employed above, (3.10), (3.11) can be 
proved and theorem 4 is validated. 

Proof of Theorem 5: First suppose that the element order condition is estab- 
lished. Obviously by the theorem hypothesis QOxcés(?) =0,¥PEP,(K) and 

on(1). $x(*;), $x(*,) are linear independent. Thus, the order of Qy is at most 
M- 3. If wEH'(K) and oy(u) satisfy the equation 


Orda (u) =0 (5.19) 


then as a consequence of the element order condition we draw out /J% w€P, (K), 
. a 

T,qgu= [TP u\ ox, Tgn= 5 1¢ tlee, Thpu= aN iT? uly, thus, du(u) =O ( 11 u) 

Then from R” =(¢,¢(u) |we/i'(K)}, we know the order of Qy is M - 3. 


Conversely, assume the order of Qy tc be M - 3, then the solution space of 
(5.19) is {¢n(P)|PEP,(K)}. Consequently, if wéH'(K), ZY we /7Y um [7¥ y=, 


Then dy (u) is a solution of (5.19). Thereupon there exists PEP,(K) such that 


oc(u) =¢,(P), therefore TP? uP, [1,gu = P\ 4, | , Tv oP , 
Os « aN ox 


ox 
and from this we can get Ju, /7? €P,(K) . This demonstrates the element 
order condition is established. Theorem 5 is proved. 


Proof of Theorem 6: We first prove (4.1). Make Jr be a set defined in 
proof a) of theorem 4. We proved that there existed a constant c(&) such 


that 
LT? wim. <e(€) twhr,+1., (5.20) 
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is uniformly established for ¥w€H"*' (K), ¥KE€J; . From reference (9) we 
know that (5.20) is equivalent to 


YEH"! (RK). | 112*wh mt <e(d) flr, +1 2 (5.21) 


uniformly established for ¥K€/;. Here w(x) = &(F7'x), +x€K. 


For ¥hEH"*'(K) the set { (13*w| KE/;} is a bounded set in H™(K). Othervise, 
forn=1, 2, ***, there exists K,¢€j. such that 

[7 2? wg >e (5.22) 
and {By,) is a bounded set, so there exists a series {By,'} and element K, of 
{By} such that By,» * By,* and thus 5, Kn’ > ; ‘Ky? G; ¥Kn’ uniformly con- 
verges to Gy. x - Consequently, in H™(K), G;, Kn' * Gy, K.» 1<j<*M. There- 
fore /7** w in H™(K) converges to nt? w. That is: ne , ‘ie bounded which is 
in contradiction to it being unbounded. Thus, { I7$*w\KEl 3 is a bounded set 
in H™(K), Moreover, for ¥A€/;, from the theorem conditions there exists a 


c(K) such that 


GE +R), UTE mie <e(K IGE 1.8 


Consequently from the resonance theorem (reference (10)) we know that (5.21) 
is established which is the same as (5.20) being established. 

Separately, using affine transform techniques we know (9) there exists a 
constant c(£) such that ¥wEH"*t'(K), ¥KE/J; we have 


af w—Pi,.4 )|whe * 
remake bryet., GOCE) Wily ey (5.23) 


From (5.20) obviously we have established for ameHnt!(K), ¥KE/J;, the 
inequality 


jo — 79 wim, Sold) pelyei., (5.24) 


Note that the conditions of the theorem have 


jo— 179 win = inf f(w—P) —79 (w—P) ha, 
PEP nth) 


qed) inf Jw Phnes.. Cold) whys 
Veli K) 
that is the inequality 


jw— 11 wim, ele) iwiner.. wel! (K) 
(5.25) 


is uniformly established forwK€ Jr. 
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For ¥KE€K., use the notation of proof b) of theorem 4. From (5.25) for 
¥oeH"*1(K), 1 =0,1,--,m, we have 


jo— TF wl, , = Oe \e— DP &\).¢ C OF |G lr+18 


ach ight) |Wieeet., =o, ir) [Wiese 


and again using (5.11) we can get (4.1). 


Similarly (4.2)~ (4.4) can be proved and using methods similar to reference 
(1) we can get (4.5). Theorem 6 is proved. 


Proof of Theorem 7: From thecrem 1, theorem 3, and theorem 4, we need only 
prove: if v,€U,(r=1,2,--) and suplo.dit.zu)<oo, then for x eeCT (FR), 

Jim Ts. se, 00) =0 (¢=0,1,2, j=1,2). 
Denote ».«<J].w,, w.€H'(Q)1Hi(Q), then using the partial integration formula 





we can get 
27,..(¢,.0.)=2 >> (@llY w, +82 [TY w,)dx 
Kex,?* 
=< {{ e(21TY w,—8,5, 19 w.—8x, 117 w,)dx 
KEK, Jf 


+{. of (Ty w,— [1 ,w,N + IT, ,e.N JN, 
+ (7? w,— 171) ,wN,— 15 ,e.N Ny jds 


+{ d.y (TY w, — 04, 19 w,)dx + [daw (lp w,— Meg.) N eds 
z 


_{ oa (119 w,—6,, 19 ~ d.9 (18 w,— lagu) N dt } 





. 


_ = {J enw, 1? w+. li_ Ni) 115 wed 


+| [1 pgt(OxpN ,-A apiN ,) ds } 
oe 


For a nonnegative integer 2, denote pi L?(K) + P (Kk) to be an orthogonal 


projection operator, P Pe L?(F) +P 2 (F) to be an orthogonal projection 


operator. Using formulae (1.4), a. 5), and (1.7) and the fact that {U. } 
passed the IPT test, we have 
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27 ,.:(¢.0%.) = 7 {| (o— Pte) (2/1 Wy —Ox, ij w.Ox, iT w,)dx 
Kek, ** 


+| @- Poe) CUTE w,— 115 ,w.N + Mia N)N, 


+ (119 w,.—11%,,.w.N,- I13,a,N,)N,jds 


” 


Tr Lex (TE WW, —Ox, iT’ we) +9Ox@ (i7? w,—Ox. hs w,) Idx 
4 


, 


+ og (OPN 1 —Ox@N 5) (12 w,—s¢@,) ds } 


- 





+ SO |) Meee (8x9 0x90 ,)ds 
Kex,’ ** 


+2 & | (p— Pte) (2N,N [Ta +(N!-Ni) 3 ,0,3ds 
KeK, FcaK?" 


+ 5 | e—Pt9)c2N.N 075, — Mm, 
Kek, 


+(Ni=N}) (113,—[15,)w,Jds 


Utilizing the Schwartz inequality, the interpolation inequality (reference 
(1)), the inequalities (3.4) ~(3.11), and the inequalities below 
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Consequently, fim Tale de) 0, Similarly, we can prove for 
¥iw0,1,2,j/1,2, limT,,(@,0,)=0 the inequality (5.26) can be gotten from 
i hs) 


(3.5) amd (3.7). (5.27) can be arrived at using methods similar to the proof 
of theorem 4. For (5.28) see reference (7). Theorem 7 is proved. 


Proof of Theorem 8: Under the conditions of the theorem, a(*,*) is uniformly 
positive. Consequently, there exists a constant c unrelated to Tt, such that 
the following inequality is established 


FTL 2, [Eu—ajes so <e{ mia 1Lu—v.leine, + sup |a( Ea, 4) —f (v9) } 
wel, hevdt, lvejct.no S (5.30) 
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From (4.5) we know that for T = 1, 2, ***, uniformly 


t 
min }£u,—vehez,20< 6 AY? \tylrios.0. (5.31) 


. se) 


For the second term of the left side of (5.30), using a method similar to the 
proof of theorem 7, we can get for +v,€U,, the inequality 


et 
jo(Eu,, v,) — f(ve) |<elodeod DO | tty Irves0 + AY! |tgir’ } (5.32) 


to be established, in which c is a constant unrelated to uy, tT. In this way 
we obtained (4.11). Theorem 8 is proved. 
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PHYSICAL SCIENCES 


INFLUENCE OF METAMORPHISM ON CHEMICAL COMPOSITION OF MACNO-FERROMICAS IN 
PRE-CAMBRIAN REGIONAL METAMORPHIC AREAS, NORTH CHINA 


Guiyang KUANGWU XUEBAO [ACTA MINERALOGICA SINICA] in Chinese No 1, 1987 
pp 27-36 


[English abstract of article by He Yixing [6320 5030 5281], et al., of 
Changchun Geological College] 


[Text] The areas studied include Anshan (Liaoning), East Hebei, Miyun (Hebei) 
and Nei Monggol. The types of host rocks of magno-ferromicas include meta- 
sedimentary rocks and meta-volcanic rocks. 


The Al,0, content and Mg/Mg/(Mg+Fe**) ratio in magno-ferromicas vary with the 
composition of the magno-ferromicas, especially the Mg/(Mg+Fe**) ratio, and 
show regular variations with the grade of metamorphien. 


The authors’ investigations show that the grade of metamorphism (mainly the 
temperature) tends to increase from Anshan through East Hebei, Miyun, and on 
to Nei Monggol. In the Shuichang mining area, East Hebei, the grade of 
metamorphism tends to increase from south to north. The chemistry of magno- 
ferromicas may provide valuable information for stratigraphic correlation, 
recognition of geological structures and ore exploration. 
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PHASE RELATIONS IN SYSTEMS Ag,S-Cu,S-PbS AND Ag,S-Cu2;S-Bi2S; AND THEIR MINERAL 
ASSEMBLAGES 


Guiyang KUANGWU XUEBAO [ACTA MINERALOCICA SINICA] in Chinese No 1, 1987 
pp 9-18 


[English abstract of article by Wu Daqing [0702 1129 3237] of the Institute of 
Geochemistry, Chinese Academy of Sciences, Guiyang] 


[Text] The phase relations in the ternary systems Ag;,S-Cu,S-PbS and Ag,S-Cu,S- 
BizS, have been studied using the silica vacuum technique. 


In the system Ag,S-Cu,S-Bi,S,; the phase relations are dominated by join-lines 
from galena to f.c.c. (Ag,Cuz-,S) and b.c.c. (CuyAgs-,S) at 500°C. With 
decreasing temperature, galena can coexist with all phases on the Ag,S-Cu,S 
join. 


There are six solid solutions and one new phase, i.e., "C" with a composition 
of Ag,.:Cuy, Bis 95,2, in the system Ag,S-Cu,S-Bi,S, at 500°C. Pavonite 
(AgBi;Ss5) contains 14 mol percent Cu,S in solid solution, but only 3.0 mole 
percent Ag,S in CuBi,Ss, solid solution. The CujBisS, s.s. and wittichenite 
(Cu;BiS,) s.s. can form join-lines with pavonite s.s., having the maximum 
content of 9.0 and 18 mole percent Ag,S. The most striking feature is the 
presence of bejaminite as a stable phase with a chemical formula of Ag,Bi,5S, 
on the Ag,S-Bi,S,; join. AgBiS, of the PbS type occupies a fairly large field, 
with a maximum of 23 mole percent Cu,S. 
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PHASE RELATIONS IN THE SYSTEMS Cu,S-PbS-Bi,S,; AND Ag,S-PbS-Bi2zS; AND THEIR 
MINERAL ASSEMBLAGES 


Guiyang KUANGWU XUEBAO [ACTA MINERALOGICA SINICA] in Chinese No 1, 1987 
pp 19-26 


[English abstract of article by Wu Daqing [0702 1129 3237] of the Institute 
of Geochemistry, Chinese Academy of Sciences, Guiyang] 


[Text] The phase diagrams of the sysiems Cu,S-PbS-Bi,S, and Ag,S-PbS-Bi,S, 
have been investigated. This paper studies the complete solid solution 

between bismuthtite and aikinite above 300°C in the system Cu,S-PbS-Bi,5,. 
The synthetic phases CuBi,S, and Cu,;BisS, have solid solution ranges in the 
ternary system with 9 and 26 mole percent PbS at the maximum, respectively. 


A complete solid solution between PbS and AgBS, divides the phase diagram 
of the system Ag,S-PbS-Bi,S, into two parts: Bi-rich and Ag-rich. All 
sulfosalt minerals and solid solutions, including pavonite s.s., lillianite 
s.8., heyovekyite and benjaminite, are on the Bi-rich side. In addition, 
bivariant relationships were found between pavonite s.s.-lillianite s.s., 
benjaminite and bismuthtite as well as between lillianite s.s.-bismuthtite 
and galenobismutite. 


Synthetic experiments using the LiC1-KCl flux technique show that when a 
minor amount of copper (less than | wt percent) is added, many Ag- and Pb- 
bismuth sulfosal minerals, e.g., vikingite (AgsPb,Bi,5559), can be 
synthesized successfully, particularly at 400°C. This was true of heyrovsky- 
ite, which has a solid solution range with 3.7 mole percent Cu,S at the 
maximum in the system Cu,S-PbS~-Bi,S,. 
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CHARACTERISTICS AND ORIGIN OF FELDSPAR MEGACRYSTS IN CENOZOIC BASALTS FROM 
SOME LOCATIONS OF EAST CHINA 


Guiyang KUANGW XUEBAO [ACTA MINERALOGICA SINICA] in Chinese No 1, 1987 
pp 37-46 


[English abstract of article by Qiu Jiaxiang [6726 1367 7534], et al., of 
Wuhan College of Geology] 


[Text] Presented in this paper are field and laboratory data on feldspar 
megacrysts in Cenozoic basalts from East China, including 76 wet chemical 
or electron microprobe analyses of feldspars, and these data are compared 
with those for 53 feldspar megacrysts from abroad. 


The mode of occurrence, physical properties, mineralogical chemistry and trace 
element distribution of feldspar megacrysts are studied using various 
approaches. 


Also presented are the authors’ viewpoints on the origin of feldspar mega- 
crysts-~a problem which has been involved in controversy for a long time 
among geologists worldwide. 
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INFRARED SPECTRAL STUDY OF COOKEITE 


Guiyang KUANGWU XUEBAO [ACTA MINERAJ.OGICA SINICA] in Chinese No 1, 1987 
pp 52-57 


[English abstract of article by Liu Gaokui [0491 7559 7608], et al., of the 
Institute of Geochemistry, Chinese Academy of Sciences, Guiyang] 


[Text] Presented in this paper are the infrared spectra of cookeite from 
Henan and Guizhou provinces, China. The frequencies, absorption intensities, 
assignments of the main bands and infrared spectra of heated samples are given. 
These spectra are different from those of cookeite reported in the literature, 
but similar to those of muscovite and illite. By comparing all these spectra, 
the spectral identifying characteristics of cookeite are proposed. The 
infrared spectrometric results indicate that the structure of cookeite may 

be complicated--there is no water in the form of "H,O”" in it and it is close 
to dioctahedral in structure, but not to tricotahedral. When cookeite is 
heated to 453°C and kept at this temperature for 30 minutes, all OH groups 
are lost and its structure begins to transform toward amorphism. When the 
temperature rises to 800°C, it is transformed into mullite, amorphous Si0, 


and Li,0. 
9717 
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DISCOVERY OF Gd-Dy-ESCHYNITE 


Guiyang KUANCGWU XUEBAO [ACTA MINERALOGICA SINICA] in Chinese No 1, 1987 
pp 66-73 


[English abstract of article by Cai Genging [5591 2704 1987] of the Institute 
of Uranium Geology, Beijing] 


[Text] Gd-Dy-eschynite has been found in the U-Th mixed mineralized coarse- 
grained biotite granite, which is associated with thorite, monazite, zircon, 
orthite, magnetite, ilmenite, bastnaesite und minor fluorite, calcite, 
apatite, etc. a:b:e = 0.3851 : 1 : 0.5381. Ten crystal forms have been 
identified, including {100}, {010}, {110}, {316}, {001}, {021}, {111}, {211}, 
{311} amd {702}. It is resin yellow, strongly resinove in luster and 
conchoidal to subconchoidal in fracture. Specific gravity = 4.716 and micro- 
hardness = 331-566 kg/cm*. It possesses brittleness and strong radioactivity. 
It is transparent in polarized light, colorless to yellow, isotropic, N > 2.20 
and grayish-white in reflected light. The difference between Gd-Dy-eschyni te 
and eschynite lies in [Y],0, > [Ce],0, and TiO, > Wb,0, in the former against 
[Ce],0, > [¥]),0, and Nb,0, > TiO, in the latter; that between Gd-Dy-eschynite 
and priorite in [Y],0,/[Ce],0, = 1.54 against [Y],0;/[Ce],0, = 5 (the former 
contains higher U and Pb than the latter); that the former contains lower Th 
but higher U than the latter; that between Ge-Dy-eschynite and sinicite in 
[¥),0, > [Ce],0, and TiO, > Nb,0, and higher U and Th in the former. Since 
Gd and Dy comprise 72 percent of the total amount of rare earth elements, this 
mineral is named Gd-Dy-eschynite. 
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STUDY OF GENETIC RELATIONSHIP BETWEEN ALABANDITE AND SPHALERITE IN DAWAN 
Zn-DEPOSIT 


Guiyang KUANGWU XUEBAO [ACTA MINERALOGICA SINICA] in Chinese No 1, 1987 
pp 74-77 


[English abstract of article by Wei Qiying [7614 4860 5391], et al., of the 
Department of Geology, Beijing University] 


[Text] Alabandite, which occurs in the sphalerite ore of the Dawan skarn-type 
zinc deposit, closely intergrows with dark, high-Mn sphalerite. Other 
associated minerals are pyrite, pyrrhotite, chalcopyrite and galena, and some 
Mu-bearing minerals, such as rhodonite, Mn-allanite and rhodochrosite, are also 
rec-ewizod. Under the microscope it is very easy to confuse alabandite with 
Sph...g:° 2. When compared with dark, high-Mn sphalerite, it has a slightly 
higher .. -ractivity and lower microhardness. Election probe and X-ray micro- 
analysis data are in good agreement with those reported in the literature. 
According to the micrographic texture, the authors consider alabandite to have 
been exsoluted from dark, high-Mn sphalerite at about 350°C. 
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CONICHALCITE DISCOVERED AT PINGGUI, GUANCXI 


Guiyang KUANGWU XUEBAO [ACTA MINERALOGICA SINICA] in Chinese No 1, 1987 
pp 84-87 


[English abstract of article by Lai Lairen [635] 0171 0088], et al., of 
Guilin Institute of Geology and Mineral Resources] 


[Text] Conichalcite occurs in the oxidized ore of Karst accumulated Au-Sn 
deposits in Guangxi. Electron probe and wet chemical analyses show (percent): 
CuO (29.45, CaO (23.32), MgO (0.22), FeO (0.08), SiO, (0.33), Al203; (0.87), 
ZnO (0.60), AsO, (41.49), H,0* (4.12), total (100.48). The formula is 

(Ca, Mg) (Cu, Fe, Zn) (As, Si, Al) 0,(0H). Ng’ = 1.799, Np’ - 1.768. Biaxil 
positive, 2V - 82°. Dnean = 4-28, Hy(so) = 190 kg/cm*. Single crystal . 
diffraction yields orthorhombic lattice dimensions a = 7.33 A, b = 9.12 A, 

c = 5.79 A. Space group P2,;2,;2,;, 2 = 4. The strongest lines of the X-ray 
pattern are 5.69(4), 4.08(4), 3.40(4), 3.12(10), 2.89(4), 2.83(10), 2.60(8), 
2.57(5), 2.066(5), 1.727(7), 1.614(8), 1.571(5), 1.469(4), 1.305(4), 
1.099(5), 1.002(4). 
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FIRST DISCOVERY OF PUMPELLYITE IN WENDUERMIAO GROUP, NEI MONGGOL 


Guiyang KUANGWU XUEBAO [ACTA MINERALOGICA SINICA] in Chinese No 1, 1987 
pp 88-91 


[English abstract of article by Xu Chuanshi [6079 0278 6108], et al., of Hebei 
College of Geology] 


[Text] Pumpellyite is found in the strata of the Wenduermiao Group, Nei 
Monggol. It occurs in Paleozoic ophiolite and is associated with glaucophane, 
lawsonite, stilpnomelane, chlorite, albite, epidote and actinolite. 


Pumpellyite is grayish green-pale brown in color, and pale yellow green-pale 
yellow brown in thin sections, showing polychromism. Ng: pale yellow brown, 
Nm: light yellow green, Np: light yellow. Ng = 1.6899, Nm = 1.6768, Np = 
1.6743, Ng - Np = 0.016, 2V(+) = 45°4. The strongest lines in the X-ray dif- 
fraction pattern are: 6.98(6), 4.68(3), 3.17(9), 2.97(10), 2.51(4), 2.12(3) 
and 1.74(2). Electron microprobe analysis shows that this pumpellyite is of a 
transition type between pumpellyite and ferropumpellyite. Julgoldite con- 
tinues to increase with increasing metamorphism. 


The discovery of pumpellyite of this type has further confirmed the existence 

of a low-temperature, high-pressure metamorphic belt in the Wenduermiao Group, 
which is of great importance for further studies of plate tectonics, petrology 
and tectonics. 
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DISCOVERY OF AQUAMARINES IN NORTHEASTERN HUNAN 


Guiyang KUANGWU XUEBAO [ACTA MINERALOGICA SINICA] in Chinese No 1, 1987 
pp 92-94 


[English abstract of article by Zheng Ruifan [6774 3843 0416] of the Institute 
of Geology, Hunan] 


[Text] Aquamarine occurs in rare metal-bearing granite-bearing granite- 
pegmatites of northeastern Hunan. It is marine blue, transparent and 
hexagonally prismatic. Its crystal ranges from 2.5 to 6 cm in diameter 

and 3 to 12 cm in length. Density: 2.71; refractive index: No = 1.5827, 

Ne = 1.5750; birefringence: 0.0077; cell parameters: a = 9.23 i, c= 9.2548 
c/a = 1.002, V = 683.27 A*. The interior of the aquamarine contains minor 
amounts of small-grained inclusions. 
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NATIONAL DEVELOPMENTS 


HANDELSBLATT VIEWS COMPUTER MERGER 
Duesseldorf HANDELSBLATT in German No 63, 3 Mar 87 p B18 


[Article by Detlef Rehn, East Asia Institute e.V., Bonn: "The U.S. Giant IBM 
Is the Model for Building Up Their Own Computer Companies") 


[Text] "China's IBM has been founded!" With these words, the Chinese press 
exuberantly greeted the merger of 67 institutions of the Beijing computer 
industry, in the middle of December 1986, to form the "Great Wall Computer 
Corporation Beijing." The new conglomerate comprises 5 universities, 4 

research institutes, service and training facilities, as well as production 
enterprises for the manufacture of computers and integrated circuits. With 
50,000 employees, capital assets of about 240 million RMB (about DM120 million)-- 
just about half the total capital assets of the computer business in China-- 

and an estimated production volume of 250 million RMB (= DM125 million), it 

is the largest computer enterprise in the country. 


Corresponding to its composition, the conglomerate covers all important 
areas of the computer industry: development and production of new hardware 
and software, training of technical people, furnishing and services, support 
for computer applications, etc. The Chinese leadership places great hope in 
the enterprise, which will take the form of a corporation. Primarily, the 
leadership assigns to the corporation a locomotive function and hopes that, 
with its aid, the still very backward computer industry of the country will 
get a start and will gradually catch up with developments in the industrial 
states, 


Reforms are indeed necessary in a branch that is decisively important for the 
Chinese modernization program. At this time the Chinese computer industry 
consists of about 130 institutions with about 100,000 employees. The centers 
are situated in the coastal areas. Besides Beijing, these centers are in 
Tianjin, Shanghai, Jiangsu, Fujian, and Guangdong. In 1985, the production 
value of the computer industry was 1.76 billion RMB (= DM855 million); this 
corresponded to only a fraction, 6.2 percent, of the total production volume 
of the electronics industry. 


For comparison: In Japan, for example, the computer industry contributes 
20 percent to the total production of electronic equipment. China can 
indeed point to a broad spectrum of its own computer developments in recent 
years, starting with supercomputers of the "Galaxy" series up to the 
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IBM-compatible “Great Wall" microcomputers. This does indeed document their 
increased technological capability, but on the other hand these positive 
developments are again put in proper perspective by other phenomena: For 
example, the Chinese computer market is almost completely dominated by 

foreign enterprises as far as large and medium systems are concerned, and even 
with microcomputers the domestic market share, according to Chinese data, is 
at most 40 percent. 


Just as serious is the fact that even where import substitution would already 
be possible, domestic computer manufacturers generally still fall back on 
foreign components. In the view of the persons in charge, this has a negative 
effect on the creation of a technological and economic base for the Chinese 
computer industry. But this means that the true local portion of value creation 
at this time makes up only a fraction of the production volume of Chinese 
computer production. 


In similar fashion, reports concerning the penetration of computers into all 
areas of the Chinese economy and society, and China's entry into the informa- 
tion age, must be regarded in differentiated fashion. By the middle of 1986, 
about 7,000 medium and small computer systems as well as 130,000 microcomputers 
were installed in China. In comparison to 1980 (2,900/600), these figures 
represent impressive rates of increase. But in this connection one must also 
see that, in a leading business center such as Shanghai, the degree to which 
computers have spread among industrial enterprises is only 10.2 percent, and 
the efficiency of utilization of the computers employed is frequently very low. 


Against this background, the Chinese leadership has developed a computer develop- 
ment strategy for the current Seventh 5-Year Plan (1986/1990). By means of this 
plan, the main weaknesses of this sector are supposed to be eliminated by the 

end of the decade, working from various sides. 


--Computer application receives top priority. In the future the relevant guide- 
line here will be that application determines the production volume. This 
approach reflects the Chinese leadership view that the competitiveness of their 
own products can be raised only with a consistent marketing direction, and that 
the desired computerization can be implemented step by step. Computers are to 
be used especially in the technical reformation of so-called traditional 
industries, for example the textile industry. Marketing direction in this 
connection also means striving for close cooperation among software houses, 
service facilities, and users, in order to secure a high use efficiency for 

the computers. 


--The product spectrum of the Chinese computer industry likewise should change 
fundamentally according to market requirements. Up to now there was a supply 

of an assortment of the widest variety of computers, which generally were 
produced only in very small numbers. By 1990, production will be more stream- 
lined. Microcomputers (also called personal computers such as e.g., the 
IBM-compatible "Great Wall" 0520) will be developed and produced preferentially; 
as need demands, minicomputers (32-bit superminis) and larger systems will also 
be produced. The objective is to reach international quality standards within 
3 years and to reach international price levels within 5 years. 
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Computer peripherals is a branch that is very neglected in China, and it too 
is supposed to gain in importance. What is planned here is primarily the 
production of technologically demanding products, which presently must largely 
still be imported: diskette drives, monochrome and color monitors, as well as 
24-pin matrix printers. 


--In view of the large technological gap of the Chinese computer industry with 
respect to the international level, China by 1990 wants to penetrate into the 
area of the fourth computer generation, which uses very large scale integrated 
circuits (VLSI). The gap would thus diminish to about 8 to 10 years. 


But implementing this project presupposes a plethora of measures. Especially 

in the area of chip development and production, the jump into mass production 

of technologically high-grade and cheap circuits and the solution of the 
associated technological and economic problems (automation of production, quality 
control, etc.) are ali required. 


--The foundation of the “Great Wall Computer Corporation" as well as further 
major mergers indicate that the Chinese leadership is particularly set on 
major conglomerates, which are supposed to cover the entire spectrum from the 
development of new systems through productionof both new hardware and software 
up to and including after-sales service. 


The responsible planners expect from these conglomerates on the one hand that 
they will overcome the system-based barriers between research and development 
and production. Another expected positive effect is that the co:.glomerates will 
be able to secure mass production on the basis of their personnel, material, and 
financial capabilities. Finally, the conglomerates should form the technological 
and economic backbone for the march into the "fourth computer generation.” 


IBM, the world’s largest computer manufacturer, is a model for the Chinese in 
founding their own computer conglomerates, especially in view of its role in 
influencing markets and products. In contrast to Western countries, the 
multiplicity of computer developments in China does not express sharp competi- 
tion and fast innovation, but rather the conceptual emptiness of their computer 
policy and the lack of businesslike market orientation. In view of this multi- 
plicity, the responsible parties hope especially that one or more of these new 
conglomerates will succeed in influencing the domestic market in such a fashion 
that its products will emerge as industrial standards. 


Quasi-private computer enterprises--which, for example, were founded by former 
scientists of the Academy of Sciences-—-represent an entirely new development in 
the Chinese computer industry and accordingly are given positive treatment. 
Because these enterprises are relatively small compared to the major 
conglomerates, and because of their operational organization, they are quite 
flexible. They can help establish entrepreneurial competitiveness. As 
spin-offs from universities and institutes of the Academy of Sciences, they 
are also useful in creating high-tech centers. 
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The import of production engineering and know-how continues to be the focus 
of the development strategy for the computer industry. But protectionist 
trends are appearing in connection with imports of finished products, 
especially microcomputers with a word length of 16 bits and less, which form 
the central focus of the Chinese computer development strategy. With these 
computers, they are striving to increase the "guochanhua," i.e., the portion 
of local production, to 80 percent of the main components by 1990. 


For West German enterprises, it follows - rom this that the best market prospects 
in the Chinese computer industry (as also in other branches) will result from 
readiness to deliver know-how and equipment. They will focus on the weak 

points of China's computer industry, i.e., drive mechanisms, hard-disk memories, 
diskettes, ail types of components, etc. The Chinese computer market is hotly 
contested. But an expected market volume of 10 billion RMS by 1990 should be 
attraction enough to enter the competition. 
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NATIONAL DEVELOPMENTS 


HELIUM PRODUCTION SITUATION DISCUSSED 


Tianjin TIANRANQI GONGYE (NATURAL GAS INDUSTRY] in Chinese Wo 4, 28 Dec 86 
pp 96-100 


(Article by Tang Wenjun [0781 2429 0193], KXinan Institute of Industrial 
Chemistry: "The Future of Helium in China as Viewed from the World Helium 
Situation") 


{Excerpts} Abstract: This paper describes the ciroumstances regarding the 
production, consumption, and supply and demand of helium outside this country, 
and also compares the levels of production and marketing in China with those 
abroad. This makes clear the gaps between the helium industry and application 
technologies in this country and those abroad. The paper proposes how we can 
learn from foreign experiences, and provides a long-range forecast and 
planning for helium production in China. 


The industry in this country that extracts helium from natural gas has 4 
certain basis, but from the point of view of scale and level of technology 
there is not only a great disparity with advanced nations, but it is also far 
from suitable for the pace of development of modern science and technology and 
of production in this country, especially in regard to its high cost and 
expensive selling price. In recent years, China has made definite advances in 
aspects of airship research, but great quantities of helium are used in doing 
so. Because the price of helium is too high, this keeps some units from going 
forward. 


The levels of production and consumption of helium in this country are very 
low. There is not only great disparity in comparison with developed 
countries, but we cannot even compare with the levels of India and Sweden. 
According to statistics, nommilitary helium consumption in developed countries 
is not less than 5 cubic km per 1 million residents (not including helium used 
for welding). Figuring from this, each 100 million in population should 
consume 500 cubic km of helium, or in this country of 1 billion, there should 
be consumption of 5 million cubic meters per year. When military needs are 
added to this, the demand for helium should be much greater than that figure. 


The price of helium in this country has always been high, the reasons being 
that the helium component in the raw material natural gas is not very high, 
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our techniques are backward, energy consumption is high, and the costs are 
high. The high price of helium has terrified many science research 
Gepart@ents, limiting the development of application technologies and is of no 
use to the godernization of science and technology. This is a probles in 
urgent need of resolution, but for the present we can only adjust the price of 
helium appropriately on the existing cost basis. The basic methods are to 
build facilities on a greater scale, use more advanced technologies, lower the 
indices of energy consumption, and greatly decrease costs. Another way would 
be to import technology, using the method of jointly funded operation with 
foreign enterprises. 


There is already a certain basis to the helius industry in this country, but 
from the perspective of scale of production, the quantities sold are just too 
much out of line. As far as technology is concerned, the cryogenic method has 
been industrialized, but because the energy consumption is rather high, costs 
cannot come down; the thin-fils penetration method is under research and 
Gevelopment, but it needs to be made more powerful to promote the 
industrialization process; the varying-voltage-adsorption technique is already 
at a certain level, but how that is to be integrated with helium extraction 
techniques awaits research and experimentation. 


Based on the actual situation in this country, I now propose the following 
recommendations and measures: 


1. Request consideration by relevant departments of the government, and in 
terms of organization, adopt measures to strengthen the management of natural 
gas resources, actively develop efforts at surveying helium resources, and 
formulate helium-protection planning that can truly be implemented. 


2. Select an area that has natural resources and build a second helius 
production base, and promote technological transformation, reductions in 
energy consumption, and reductions in cost to greatly reduce the selling 
price. 


3. Strengthen capacities and resolutely pay attention to the research and 
development of new technologies, especially the membrane separation 
technology, the (Kongfen) helium extraction technique, etc. 


4. Establish a Chinese industrial gas company, under which set up a helium 
branch company to uniformly manage all industrial gases in this country and 
uniformly manage helium resources, production, and marketing, which will 
vigorously promote the development of the helium industry in this country. 


In China, we are currently working for the four modernizations, and as the 
wave of the new technological revolution is coming to engulf us, this 
strategic material that is helium has its own particular importance. The 
dominant position of American helium can be maintained for 30 to 50 years, and 
the helium production capacity of the Soviet Union has now vaulted into second 
place in the world. When to that is added the natural resources of Poland, 
the countries of Eastern Europe can satisfy the current and future demands for 
helium. China's economic leap is still in its initial stages, and based on 
the experiences of countries like Japan, future demand for helium is certain 
to grow 100-fold. Looking at this from the perspective of the situation in 
this country regarding natural resources and technological development, if we 
do not now pay close attention to this problem, it is certain to bring serious 


consequences in the future. 
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NATIONAL DEVELOPMENTS 


CONFERENCE LECTURE ENCOURAGES WORK IN SOFT SCIENCES 


Tianjin KEXUEXUE YU KEXUE JISHU GUANLI [SCIENCE OF SCIENCE AND MANAGEMENT OF 
S&T] in Chinese Wo 12 Dec 86 p 1 


{Article prepared from lecture by He Zhongxiu [01459 6988 4423): "I Hope that 
the Cause of the Soft Sciences Will Flourish and Develop in China") 


(Text) Today, sore than 300 workers in the soft sciences have come from ali 
over this country and from sore than 10 nationalities to gather here to 
discuss the problems of youth in soft science research. This occasion gakes us 
all tremendously excited. 


When I say that "this makes us tremendously excited,” this is no exaggeration, 
and everyone feels the same. More than 2 months ago, the majority of those 
participating in the All-China Soft Sciences Research Working Conference were 
middle aged or older. More than 80 comrades spoke and many achievements were 
exhibited, all of which unequivocally proved that soft science research truly 
has a special and important role to play in the socialist sodernization of 
this country. A lecture by Vice Premier Wan Li elicited a great response 
throughout this country. At present, everyone is calling for the creation of 
conditions in which a new generation of workers in the soft sciences can 
mature. As we open our conference here today, more than half of you are young 
comrades in your twenties and thirties, and how can this not excite us? 


Someone asked, what is the aim of calling this conference? It is just in hope 
that we can create the conditions for our young comrades that will provide the 
opportunity of a gathering in which to exchange ideas, so that we sight 
promote the development of the cause of the soft sciences in this country. 


There are a minority of middle-aged and elder co-workers among us, and we came 
to the practice of the soft sciences after the breakup of the "gang of four." 
We have labored for nearly 10 years, and gost of us are in our forties and 
fifties. I am sure you all remember that at that time any proposal to engage 
in the soft sciences was not understood: what are "soft sciences"? That is 
too obscure and cannot solve any real problems; when you think about it, 
distant water cannot quench a nearby thirst; and it was even believed that 
that was "not a real occupation.” But it was in just that kind of atmosphere 
that some comrades said: "Let us proceed from the actual situation in this 
country. Let us study the role of science and technology in the economy and in 
society. Let us explore the rules for the development of science and 





technology to understand the basis for development principles and governmental 
policies regarding efforts in science and technology. Let us review the 
fundamental experiences of successful management, and let us ascend to the 
heights of scientific principles to guide our practice. This will allow 
science and technology to serve the development of socialist sodernization, 
which is, after all, the pressing mtter." They said that “even though they 
are soft sciences, they will still require hard effort," and “even if kept 
waiting for 10 years, I will still carry on with China's soft sciences!" 
Classes have gone beyond renouncing the pen for the sword and have instead 
"renounced hard for soft," “left specialties to take up management," and have 
resolutely taken up a life work in the soft sciences. 


Today, research in the soft sciences is not likely to be seen as "not a real 
jod.* We at this conference have the support of so many leading units, the 
@ass media, and research organizations that these are then the conditions 
created for young skilled personnel in the soft sciences to make their aark. 
But we cannot say that for these reasons there will be no difficulties for the 
young in studying and researching the soft sciences. There are in fact 
various new obstructions, as for example, “Are the soft sciences things that 
you young people can do?,”" and “Getting involved in fads!" ... are often-heard 
opposing views. When some comrades received the invitations to this 
conference, their unit would not support them, so they could not come. I have 
heard that some of the younger comrades at the conference today have come on 
their own expenses by requesting a rotation of their holidays. They have made 
a long, arduous journey to seek teachers and visit with friends, and they have 
not flinched in the face of adversity, which deeply soves us. 


The soft sciences are necessary for our modernizing development. Comrades who 
are determined to take up the cause of the soft sciences gust have the spirit 
to be persistent and dauntless and to not be swayed by repeated setbacks. I 
believe that both we middle-aged and elder comrades and our young comrades 
should have this drive: be resolute, fear no pressure, and persevere; ve 
tenacious, be not in awe of any hardship or danger, and you will be able to 
persist; be enthusiastic, take no advantage of anyone, and be willing to 
commit yourself; be true, engage in no fraud, and seek real results. As long 
as the direction is correct and we have these four determinations, we can be 
all-conquering. 


In recent years, there have appeared in all fronts throughout the country a 
group of youth who have set their minds on soft science research, and who are 
full of life and vitality. They have their own ideas and their own 
determination, and have produced many achievements, but there has been little 
in the way of lateral relations, and they have yet to constitute a contingent. 
If we oan create the conditions for thes, can provide an opportunity of a 
gathering in which to share ideas, then not only will we allow them to become 
known to each other, to strengthen their relations, and to form into a 
contingent, but also we can allow those who are middle aged to meet the next 
crop and have each gain from the other. With mutual cooperation, the cause of 
the soft sciences is certain to flourish even more. 
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NATIONAL DEVELOPMENTS 


CHONGQING UNIVERSITY RAISES OW! RESEARCH FUNDS 
Beijing GUANGMING RIBAO in Chinese 20 Jan 87 p 1 


[Article by reporter Li Jiajie [2621 1367 2638]: "“Creases Competitive 
Environment, Encourages Faculty to Go Off-Campus to Open Up Sources of 
Funds for Scientific Research; Lest Year Chongqing University Raises 

Over 13 Million Scientific Research Punds On Its Qwn; Exceeds Overall All- 
University Scientific Research Funds by 90 Percent™] 


[Text] A fundamental change has taken place in Chongqing University's 
source of scientific research funds. Statistics up to the end of 1986 
indicate that the self-raised funds exceeded overall all-university 
scientific research funds by 90 percent. 


An important reason why the university was able to increase self-raised 
scientific research funds by such a large amount is that the school did 
its utmost to create a competitive environment for the faculty and when 
allocating scientific research funds provided by the State it was not 
unbiased. Whichever department obtained more results received more 
scientific research iunds next year. 


Another important reason is that the schocl was fully aware that the 
faculty are in the position of masters of the school. As long as the 
faculty firmly believe that they can undertake a scientific research topic, 
even if it will be a big risk, they will firmly support the faculty 
setting up the topic. After obtaining results for such risky topics, 

they also firmly support the faculty going into the factories for 

carrying on further applications research. Through just this one move the 
university approved investing 700 thousand yuan in resources to mount 

50 topics and won back from the factories 240 topics and 6 million yuan 

in scientific research funds. 


In those instances when it cannot be commercialized, research results 
which win prizes off-campus, the school always awards the faculty school 
prizes and encourages the faculty to go off-campus to find sources of 
funds for scientific research. 


In 1986, the state appropriated only | million yuan in scientific research 
funds to Chongqing University but the faculty themselves found over 13.5 million 
yuan. 
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NATIONAL DEVELOPMENTS 


BEIJING JOINT RESEARCH-PRODUCTION ORCANIZATIONS PROLIFERATE 
Beijing GUANCMING RIBAO in Chinese 20 Jan 87 p 2 


[Article: "Seven Types of Joint Research and Production Appear in 
Beijing Area; Joint Research and Production Organizations Expand to 
450 Including 181 Which Have Spread Beyond the Region”) 


[Text] According to information from persons concerned, joint scientific 
research and production in Beijing Municipality is showing new develcpment 
momentum. The number of joint organizations is constantly increasing 

and the area covered is becoming larger. In addition, there is a tendency 
to develop from casual forms to semi-close and very close forms. Up to 
the present joint research and production organizations have increased to 
450, 181 of which extend across regional boundaries. 


Concerned persons also told reporters that 7 models of joint research and 
production have appeared in Beijing Municipality. 


Long-term Comprehensive Technical Cooperation Type 

As market competition increases, to improve their technological development 
~ lity and absorb into their own sphere social and science and technology 
forces in a variety of forms, many enterprises have formed lonyg-ters 
cuaprehensive technical cooperation relationships with scientific research 
units and institutions of higher education. This type of cooperative 
relationship generally involves the enterprise supplying the topic and 
certain conditions and the scientific research unit is responsible for 
attecking the problem and supplying results. The enterprise arranges for 
intermediate testing and putting it into production. This method of close 
linking of research, intermediate testing, and production has shortened the 
research cycle so that scientific and technological results are rapidly 
converted into productive forces. For example, the polyacrylic high effect 
catalyst successfully developed by the Yanshan Petrochemical Company in a 
three-way cooperation with the Academy of Science's Institute of Chemistry and 
Institute of Chemical Engineering in just a year passed (Xiaoshi) appraisal 
and in terms of quality achieved international advanced levels. They are 
now constructing an intermediate testing shop at the Yanshan Petrochemical 
Company and after going into production it can save 1.8 million dollars 

a year. 
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Technical or Intellectual Shares Type 

Treating their scientific and technical results or intellectual work as 
investment in an enterprise, scientific research institutes have increased 
earnings and rationally divided it among the two parties. This is a new 
form which has developed on the foundation of past technical transfer. It 
avoids the difficulty of the enterprise having to pay transfer fees before 
it has earned any income and at the same time it opens for the scientific 
research institute a rather stabie source of funds and is beneficial for 
mobilizing the initiative of both sides. 


Joint Venture Type 

This consists of a plant managed jointly by a research institute and a 
relevant enterprise, the leadership and management is by agreement between 
the two sides and both sides provide personnel. The research institute 
provides the technical results and some of the technicians, and some 

are also responsible for marketing and market development; the enterprise 
supplies the plant building, equipment and labor, and the earnings are 
divided by agreement. The experience of many enterprises proves that once 

a joint venture relationship is established with a research institute it's 
like adding wings to a tiger. For example, the Haidian Agricultural Machinery 
Manufacturing Plant had been losing money for several years and was on the 
verge of closing down. After joining in a venture with the Iron and Steel 
institute they began producing a special type of alloy steel tubing it 

was in great demand. In 6 years it created a value of production of 

23 million yuan, last year the value of production grew 14-fold over 

the period before the joint venture and profits grew 11./7-fold. The profits 
handed over to the state increased 45-fold, and labor productivity 

increased 13-fold. 


Complete Technological Development and Engineering Contracts Type 

By focusing on development of complete technology and engineering contracts, 
this type combines research, design, manufacturing, and construction units 
to form a new technological economic entity which has legal standing. 
Although there are not many of this type of joint organization at present, 
ihey are rather special and have very good prospects for development. 

For example, the China Recreational Vehicle Joint Vompany in which the 
Fersous Metals Design Institute took the lead and in which over 20 enterprises 
in Beijing and other areas participate, since its construction it has 
provided 28 provinces and municipalities with over 1,300 recreational 
vehicles in more than 80 types and has saved a great deal of foreign exchange. 


Industrial Comprehensive Development Type 

Comprehensive development focussed cn one industry combines laterally 

research and production units which are scattered in other systems. Generally 
a casual type alliance is formed first then as needed many semi-close and 
close type alliances are formed. For example, the dried and fresh fruit 
improved variety breeding technology and deep processing connected alliance 
formed by 12 units, including the Academy of Agricultural Sciences’ Orchard 
Institute, Forestry Bureau, Agricultural Academy and some counties and 
districts has formed 13 joint organizations from formulating dried and 
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and fresh fruit comprehensive development plans and beginning to promote 
scientific and technical results, has established municipal, county and 

key rural specialized fruit tree centers, has propagated 5,100 m of 
nursery, 25,000 mu of demonstration orchards, and a demonstration system of 
20 million jin of storage and processing capacity. Direct income from just 
the excellent nursery stock provided in 1986 reached 2.6 million yuan. It 
played an important role in promoting fruit products production in the 
Beijing suburbs. 


Jointly Run Collective Scientific and Technical Development Agency Type 

In the past few years a large group of scientific and technological 
development agencies owned by collectives and individuals have appeared 

in the Beijing area, mostly in the Haidian district which has over 100 of 
them. Many of them are run jointly by the Haidian district, research units 
of the Chinese Academy of Sciences and various departments of the State 
Council, and institutions of higher education, and there are also many 

which are operated by scientific and technical personnel who have voluntarily 
thrown away “iron rice-bowls" and formed partnerships. Sitong, Kehai, 
Haihua, and Xintong have already obtained interesting results which have won 
praise from various areas of society. The value of output of these 

few companies reached an unprecedented 200 million yuan in 1986. 


Joint Scientific Research and Production Group Type 

This is a new type which has suddenly appeared. It generally relies on a 
tap industry and has as -ts goal developing new technology, new products, 
new systems, starting new industries, and bringing along old industries, 

it breaks through boundaries, develops regional advantages, and links 
research, production and applications units into a group company. For example, 
the Great Wall Computer Group Company recently approved by the State Council 
and has 67 participating units is new type enterprise group which is a 
combination of research, development, production management, applications 
service, and training which encompasses 50,000 employees and over 1,500 
scientists and technicians. Its establishment will promote development of 

a northern electronics industry base and provide new technological material 
foundation for development of the economy of the capital and even entire 
nation. 
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NATIONAL DEVELOPMENTS 


TECHNOLOGY MARKET SURVEY 
Tianjin JISHU SHICHANG BAO in Chinese 2 Jan 87 p 1 


[Article: "Like the Constancy of the Moon, Like the Rising of the Sun; 
A Technology Market Survey Report"; first paragraph is source-supplied 
introduction] 


[Text] Editor's note: In the last half of last year, the National 

Technology Market Survey Directorate organized a survey group jointly with 
the State Science and Technology Committee, the State Council's Legislation 
Bureau, State Industry and Commerce Administrative Mayagement Bureay, and the 
Tianjin Municipal Science and Technology Committee carried out a nationwide 
technology market survey. They visited 4 provinces, | district, and 12 cities 
(Liaoning, Heilongjiang, Zhejiang, Hunan, Guangxi, Shenyang, Liaoyang, 

Harbin, Tianjin, Shanghai, Hangzhou, Guilin, Nanning, Zhuzhou, Liling, Yiyang, 
Changsha) and in 4 months visited research institutes, plants and permanent 
technology markets and held dozens of informal discussions of different 

types. Their survey is offered to the reader below. 





Quick Action on the Reform cf the Science and Technology System Has Had 
a Big Impact; Science and Technology Market's Have Put Down Firm Roots in 
China 


Since the CPC Central Committee announced their decision on reform of the 
science and technology system, the commercialization of technological results 
and opening up of the technology markets has moved rapidly nationwide. 
According to incomplete statistics from 34 ministry committees, 29 provincial, 
municipal and autonomous district and 7 cities with province-level ecoromic 
decision-making authority, there are over 5,000 agencies engaged in tchnology 
trade and the volume of trade actually conducted in 1985 in the technology 
market was over 2.3 billion yuan. In the reform of the science and technology 
system, the technology market moved the fastest and was influenced the 
greatest. It has put down firm roots in China. 


However, the technology market is something new. Since the last half of 
1985 it has encountered conflict several times and since they did not 
fully understand the documents from upper echelons, some departments and 
local areas were not clear about the dividing lines between some policies 
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and confused opening up technology markets with the incorrect style of 

party and government agency engaging in trade and treated them as "fly-by- 
night companies." Some scientific and technical development and technological 
product management agencies were subjected to various kinds of pressure. 

The Central Committee promptly discovered and corrected these phenomena so 
that the vigor of nationwide technology markets would develop again 

and the technology market is now steadily advancing. 


Currently, nationwide various types of technology trade exchange agencies 

have been springing up like bamboo shoots after a spring rain. Science 
committees, economic committees science associations, industrial societies, 
industrial offices, academies of sciences, education committees and democratic 
parties, mass organizations, and civilian organizations have established their 
own technology trading agencies and various areas have formed preliminary 
multi-level, multi-channel, multi-form technology trading network. Various 
“armies” have formed their own powerful “front armies," which have both 
horizontal alliances and mutual competition and promoted the development 

of technology markets. 


According to Tianjin's statistics of a year ago, the entire municipality 
had over 570 various types of technology grading agencies, involving over 
60,000 persons, distributed in various industries and trades and cities, 
districts and counties to form a close alliance of science and technology 
and the economy and an important force for promoting economic development. 
The technology development business agencies and branches now established 
in Shenyang number over 6,000. The Shanghai technology market is 

nicknamed "8th Route Army": i.e., Shanghai Technology Development Exchange 
Center, Shanghai Industrial Foundation, Starfish Technology Development 
Company, Science and Technology Consulting Service Center, Shanghai 

Higher School Science and Technology Service Center, Employees Technology 
Association, and Civilian-run Technology Development Agency. With their 
own solid strength they form an enormous rank. According to incomplete 
statistics from markets in Guangxi, over 200 technical trade agencies 

for technical development and technical consulting have been established. 
Hunan Province has also established 574 technical development, technical 
service, and technical consulting agencies. The provinces have preliminary 
technical market networks of different levels. 


Technology trade forms are now developing profoundly in multiple functions 
and multiple directions. 1) Trmsferring pure technology to lower levels 

as a package. For example, the Tianjin Technology Association and a light 
industry bureau helped the Dagang Tool Works No 6 to start an imported loop 
cable production line, and accepted and completed such mission as translating 
blueprints, checking equipment, installation and testing, and connecting 

the machinery and producing and ensuring that the production line would 

go into production and operate normally. 2) Intermediate service agencies 
developing in-depth service, gradually changing the method of purely picking 
up service expenses. The Shanghai Science and Technology Development Exchange 
Center stressed intensifying the creation of “talent banks" and management 


95 








of civilian-run science and technology agencies, exchanged various types 

of talent so that civilian-run agencies could make a contribution to 

national economic construction. 3) Vigorous development of lateral 
connections between scientific research and production. Through technology 
market mechanisms many areas have established allied cooperative organizations 
of scientific research, education, and production in many forms and according 
to incomplete statistics, over 9,800 allied organizations of scientific 
research and production have been established nationwide. 4) Developing 
technological cooperation internationally. Some have developed series of 
products through joint investment with foreign factories and merchants, some 
have arranged technology exports, some units in cooperation with foreign 
merchants have entered the international market with advanced technology. 

5) Assisting enterprises to develop technological demonstrations so that 

they will change from reckless importation of technology will change to 
consultation and verification and scientific policy determination which 
promoted enterprise development. 


Stimulate the Economy, Science and Technology Relieve Inadequacy, Train 
Talent, the Social Benefits Produced by the Technology Market Are Clear 


Practice of the past few years has shown that development of the technology 
market has promoted technology transfer, stimulated the local economy, 

and brought obvious social benefits and economic benefits. On the basis of 
the Sharghai Municipal Science and Technology Development Exchange Center's 
tracking analysis of 104 technology contracts signed in 1985, each yuan 

ef contract increased the value of production by 16 yuan, increasing 

income 8.6 yuan, scientific and technical personnel receive 1 yuan, 

and the state increased revenue by 36 yuan. Seen from the perspective of the 
31 Center projects checked and accepted in 1985, invested expenses were 

7.5 million yuan, value of annual output was 62.792 million yuan, 
equivalent to 2.6-fold the total investment. The results of some single 
projects the benefits were even higher. The instances of factories and 
mining enterprises in other areas such as Liling City and the Viyang 
District in Hunan and even Huayuan County in western Hunan making up 
deficits and increasing surpluses, coming back to life from the dead, 
increasing outpuc and income and escaping poverty and becoming rich 

through importing technology and development new products are too numerous 
to mention. 


News also keeps coming in of the victory of the technology markets in 
expanding science and technology to relieve poverty and to organize and 
implement “spark plans.” For example, around a “spark plan", the Hunan 
Provincial Science Association's technology market broadcast technology 
supply and demand information, held 2 large information meetings, printed 
4 volumes of material, issued over 8,000 items of information and 
abstracted 417 information items on “short even and quick” technology 
results. The Yiyang District technology market actively supports the 
Yiyang County Microbiology Institute's implementation of the state's “spark 
plan" new technology for cultivating black tree ears, it has improved 
technical training, set up demonstration plants, and has helped 
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10 cultivation specialists to escape poverty and become rich. In the past 
few years, the Zhuzhou City Technical Association has sent out 

118 specialist technical service brigades for technical assistance to over 
80 enterprises in Guizhou, Yunnan, Guangxi and the border areas of western 
Hunan, they saved 31 enterprises which were on the verge of closing, 

with new increases in the value of output of over 70 million yuan. 

The Guangxi Science and Technology Development Center undertook such 

tasks as the formulation and implementation of a “spark plan” for the 
entire district and the technical development of the mountain region. 

They helped the counties in the mountain region formulate and develop overall 
plans, arranging technical development projects for 17 mountain regions 
with an investment of 16.78 million yuan. 


To suit the needs of technical market development, technical markets have 
trained a large number of technical administration management personnel 

ind various types of specialist technicians. The Shenyang Technical Market 
Research Association has been assigned by the municipal science committee 

the task of standardizing and giving specialized training to middle level 
personnel in various agencies and technical market management personnel city- 
wide, beginning in June, 1986, and has already held 4 classes and trained 

240 persons, and by the end of the year will have trained over 1,000 persons. 
In the past few years, the Tianjin Municipal Science and Technology 
Consulting Service Company has trained over 2,000 specialized technicians. 
The Yiyang Regional Technical Market in Hunan has held 36 specialized 
training classes and trained 2,670 technicians. The Zhuzhou City Employees 
Technology Association has begun specialized technical training in which 
over 10,000 persons have participated and over 600 lectures on technology 
have been held providing a solid foundation for technical market development. 
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NATIONAL DEVELOPMENTS 


FATIGUE STRENGTH TEST ROOM DEVELOPED 


~ 


Beijing KEJI RIBAO in Chinese 16 Jan 87 p 2 


[Article by Wang Wergvo [3769 5898 0948]: “Aircraft Fatigue Strength 
Test Room Constructed; Brings Glad Tidings For Aviation Industry and 
Air Units"] 


[Text] China's first aircraft fatigue strength test room of 1980's level 
was formally established on 30 December 1986 at an Air Force Research 
Institute. 





The aircraft fatigue strength test room is a new facility for testing and 
studying the life and reliability of air force aircraft currently in 
service. In the past, the “lead” aircraft method was used to determine 

the useful life of domestically manufactured aircraft, i.e., first an 
aircraft was test flown and if no trouble occured with the aircraft in 

300 hours then the other aircraft of the same type were provisionally 
considered to have a useful life of 300 hours. Then the “lead in” aircraft 
continued to be flown. To a considerable degree this method was a risk. 


The birth of the aircrcft fatigue strength test room is glad tidings for 
the aviation industry and flying units. In this test room various types 
of aircraft can be fastened to the test frame completely and various 
inflight load maneuvers simulated. Deformations and cracks produced 

in a part measured by monitoring of the aircraft through 512 channels 
can be accurately recorded in the control system display and the useful 
life of the aircraft determined. 


The test room is not only suited for use in measurement and research 

on various types of aircraft but also has broad civilian use for testing 
force-bearing parts of motorized vehicles, metallurgical machinery, 
agricultural machinery and small ships. 
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NATIONAL DEVELOPMENTS 


CAS SHANGHAI EARNS FOREIGN EXCHANGE THROUGH EXPORTS 
Shanghai WEN HUI BAO in Chinese 11 Jan 87 p 1 


[Article by Zhou Yuan [0719 0955]: “Open up the International Technology 
Market, Establish Joint-Venture Companies; Chinese Academy of Sciences Branch 
Creates New Path Using Science And Technology To Create Foreign Exchange; 
High Tech Products Marketed All Over The World, Last Year Export Volume 
Exceeded Import Volume" ] 


[Text] The research institutes of the Shanghai Branch of the Chinese Academy 
of Sciences have opened up the internationai technological market and have 
signed with foreign plants and merchants for cooperative development projects 
and established joint venture companies. Many high tech products are marketed 
throughout the world, last year creating a record of export volume exceeding 
import volume. This new road using science and technology to create 

foreign exchange clearly has a bright future. 


In some cooperative research and development project, the achievements of 
the Shanghai Institute of Pharmacology are striking. This institute has 
signed 7 contracts with 6 pharmaceutical companies in Japan and the United 
States to develop new medicines cooperatively. They include antineoplastic 
drugs, gastric ulcer drugs, new antibiotics, and health beverages that can 
treat melituria. Generally, after the other side has screened the preliminary 
research results provided by our side, they have priority production rights 
and apply for a patent in the relevant country but the inventor's rights 
belong to us. Foreign pharmaceutical companies annually provide scientific 
research expenses to the Shanghai Institute of Pharmacology, and once patent 
results have been spread, the profit is divided between the two sides. 
Clinical tests on an anti-cancer drug jointly developed by the Institute 

of Pharmacology and the Japanese Zenyaku Co., Ltd joint stock company have 
already been completed and 10 companies in Europe and America have applied 
for patents and many famous foreign pharmaceutical companies have expressed 
interests in buying it. g§o that Chinese medicines can better benefit 
mankind, this institute and the Japanese Yamanouchi joint stock company has 
sigu d a joint 5-year development contract and has also made advances. 


Since last year, in Shanghai the first joint ventures companies between 


institutes and foreign merchants have been born in the Chinese Academy of 
Sciences system. Last July, the Shanghai Institute of Technical Physics 
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set a precedent in establishing a joint venture with a Japanese company-- 
(Nisaila) Company Ltd.--in which both sides will provide technology, 

equipment and funds, and most of the products will be sold in the international 
technology market. This joint venture is curretlyinstalling equipment 

in Shanghai and is about to go into trial production, The Shanghai Silicate 
Institute has established a joint venture with the Yongqing Company of 
Singapore primarily to produce new types of energy-saving insulating 

materials and will construct a petrochemical industry enterprise in 

Southeast Asia to take care of pipeline design and materials supply. 


From the recent data of the Shanghai Branch of the Dongfang Scientific 
Instruments Import-Export Company of the Chiense Academy of Sciences, the 
Shanghai branch has thoroughly turned around the past situation of “only 

use foreign exchange, not create foreign exchange". Last year's export 
volume exceeded import volume by one-third. The BGO crystals of the Shanghai 
Silicate Institute are the highest in output volume in the world and 

a high tech product of the highest quality. The large order from the 
cooperative group of the Western European Nuclear Center headed by 

Professor Ding Zhaozhong [0002 5128 0022] will be used in key parts of 

the world's largest electron-positron collider. Last year the Silicate 
Institute produced over 2,000, creating over 2 million yuan in foreign 
exchange. The Shanghai Optics and Fine Mechanics Institute's YAC crystal began 
to enter the U.S. market and is used by world famous scientific instrument 


companies. 
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